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ABSTRACT 


The purpose of this research was to investigate the effect of tidal- 
period internal wave vertical motion on oceanic mixed layer dynamics, 
and to discern the effect upon the diurnal evolution of thermoclines, 
as observed during MILE (mixed layer experiment), Vertical advection 
was added to a one-dimensional bulk model of the mixed layer with an 
assumed linear in z and sinusoidal in time dependence. The rate of 
mixed layer deepening was therefore due to the combination of vertical 
mction and entrainment. The first significant result was the finding 
that the interaction between vertical mixing and vertical motion 
depended upon the wave frequency and its phase relation to the diurnal 
heating cycle. Second, linear and non-linear interactions of the wave 
induced vertical motion with the cyclical boundary conditions can 
generate two-dimensional (x-z) structure in the near-surface temperature 
field of an initially horizontally homogeneous ocean under the influence 
of horizontally homogeneous surface boundary conditions. Finally, this 
advective interaction increases the utility of the mixed layer model in 


Single station forecasting, 
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T. INTRODUCTION 


In recent years much study of internal waves has focused on the 
linear and non-linear wave interactions. Models for wave prediction 
have been developed and tested to match the theory. There have been 
attempts to measure the internal wave fields in situ. Some of the 
most useful of these have been the observations of temperature and 
pressure field variations with depth, providing a basis for computa- 
tions of the vertical motions and interactions involved. These studies 
and analyses of data demonstrate the importance of the internal wave 
field in upper ocean and near surface dynamics. 

Roth et al. (1981) presented clear evidence that internal waves 
are not as critical in an analysis of the dynamics of mixing at depth 
because they have little variability. However, near the surface, 
where they interact with other phenomena such as cyclical heat flux, 
wind stress, and other boundary conditions, their effect can be 
considerable when dealing with the budgets of heat and momentum. 
Holloway (1980) pointed out that internal waves are not weak waves 
(in the sense of weak resonant interaction theory) and have a high 
impact on mixed layer dynamics down to and including turbulent scales 
of motion. Pinkel (1981) showed that the near-surface wave field has 
a profound effect on the vertical velocity field in the upper ocean. 

Some mixed layer experiments such as those by Linden (1975) and 
Kantha et al. (1977) have shown that internal wave radiation may play 


an important role in controlling the rate of deepening of a turbulent 


1] 





mixed layer in a stratified fluid. Displacements induced by the turbu- 
lence in the base of the mixed layer served to form internal waves 
which removed kinetic energy of motion from the mixed layer. Since 
this energy is also important in the entrainment process, these induced 
internal waves, of relatively high frequency, serve to retard deepening 
of a mixed layer. Townsend (1968) has derived an expression for the 
rate of generation of the internal waves and their energy by a steady 
disturbance field which is advected relative to a stratified layer at 
a constant speed. The theory was extended by Bell (1978) to cover 
general space-time dependence of the disturbance field. These works 
indicate that internal waves could be a means of retarding the growth 
of mixed layers, through the removal of the kinetic energy in the layers. 
Such studies all indicate that the presence of internal waves should 
be taken into account as a possibly important factor in the mechanical 
energy budget which controls mixing. It is hypothesized that internal 
waves which have periods comparable to the dominant time scales of the 
atmospheric forcing (surface heat flux and wind stress) will be the 
most likely to result in an interaction with the entrainment mixing 


process in terms of a positive feedback (increased mixing) effect. 


A. BACKGROUND DATA ANALYSIS 

During the period August-September 1977, a mixed layer experiment 
(MILE) was conducted at Ocean Station P (55° N, 145° W) in the North 
Pacific. The depth of the mixed layer varied between about five and 
thirty meters during the period of the experiment. Temperature was 
continuously observed at fixed depths for a period of three weeks. A 
spectrum analysis of this data (See Figure 10 in later section) shows 


some important features. 
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First, a major spectral peak was at about 12.79 t .02 nours, with 
secondary peaks at 6.62 * .13 and 3.94 * .06 hours. Lower order 
spectral peaks are at 4.9 hours, and at several frequencies having 
periods under 3.5 hours. The peak at 12.79 hours is clearly asso- 
ciated with the semi-diurnal tide (12.45 hrs.), with the other two 
peaks at 6.62 and 3.94 hours probably being higher harmonics. At all 
three depths of the analysis it is clear that the semi-diurnal internal 
tidal wave, with its harmonics, is the dominant cause of vertical 
motion at the base of the mixed layer. 

Second, the tidal peaks, including the harmonics, are very coherent 
vertically throughout the depths analyzed. They are present in all 
three depths analyzed, below the layer, and in the mixed layer, with 
very little frequency or phase shift. 

Finally, the spectral peaks at other than tidal frequencies are 
not vertically coherent. Also, there is a large variation in non-tidal 
spectral densities for observations above and below the thirty meter 
depth--the top of the seasonal thermocline at this time of year. This 
Suggests that something has happened to the higher frequency waves 
near the surface making the interpretation of their effect on mixed 
layer dynamics hard to judge from simple theories, However, at the 
same time it suggests that the effect of tidal internal waves may be 
quite easily discerned, 

These results confirm the work of Pinkel (1981) whose studies 
showed that of the wave fields in the upper ocean, the tide and its 
harmonics have the greatest vertical coherence. 

Studies by Simpson and Paulson (1979) demonstrated the variability 
of the temperature structure in the upper ocean during the POLE experiment. 
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Shook (1980) showed that these variations. were a response to vertical 
and horizontal advective processes. In a study of Pacific Ocean ship 
track temperature data, Butler (1981) found definite cyclical varia- 
tions in the sea surface temperature structure, 

From these studies of actual data, it is obvious that (1) the tidal 
period internal waves and their harmonics are an important source of 
energy in the upper ocean, and (2) surface layer temperature patch- 
iness exists in the ocean, and may in part be generated by interaction 


of mixing and internal wave motion, 


B, ADVECTION IN HODELLING 

To date little has been published concerning the interaction of 
tidal-period internal waves and entrainment by the mixed layer. Je 
Szoeke (1980) has discussed advection in the mixed layer as a result 
of horizontal variability of the surface wind stress. He showed that 
advection, upwelling and downwelling, was an important parameter in 
mixed layer dynamics, resulting from the boundary forcing of the wind 
stress. Cushman-Roisin (1981), Garwood et al. (1981), and Price (1981) 
also demonstrated that advection can play an important role in the 
dynamics of the mixed layer, but the interaction with internal gravity 
waves is not treated. De Szoeke and Rhines (1981) discussed coastal 
upwelling at length, but the advective motion was conceived as a 
necessary factor to achieve a realistic balance for surface heating, 
transport, and mixing, in order to produce observed coastal temperature 
fields. Price et al. (1978) included vertical advection in the test 
of a mixed layer formulation during storm conditions. They found that 
the phase dependence between wind stress and mixed layer deepening, 
and the mixed layer deepening itself were reasonably simulated. 
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Within this discussion of a limited amount of relevant literature 
on the advective process it has been shown that vertical advection 
must play an important role in storm response, coastal circulation, 
and fronts, No one has examined the role of internal wave motion, on 
the tidal scale, in the open ocean where the geometry and boundary 
conditions are rather one-dimensional. 

In possibly the most relevant previous effort, Spigel (1980) 
looked at the coupling of internal wave motion and entrainment in a 
two-layer lake model, having a mixed layer which is based upon the 
one-dimensional turbulent kinetic energy budget. Time scales of 
internal wave motion, wave decay, and entrainment were widely 
separated, Using a seiche effect in a basin, he demonstrated that 
internal waves did in fact play an important role in the entrainment 
zone. 

With this background it was believed that the time has come for 
an investigation of the effects of internal vertical motion in the 
oceanic mixed layer. The motion is there. Observations snow that 
it dominates Eulerian measurements in the upper thermocline. The 
importance of this motion to mixed layer evolution needs to be 


established, 


Ce SCOPE OF THIS STUDY 

The scope of this initial study shall be confined to a limited 
number of wave frequencies. There is insufficient time to concentrate 
on all the types of internal vertical motion, and all of the possible 
couplings and interactions in this thesis. The purpose of this initial 
research was to investigate the effect of dominant tidal-period internal 
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wave vertical motion on oceanic mixed layer dynamics, and to discern 
the effect upon the transient thermocline as observed during the mixed 
layer experiment. The effects of the vertical motion are explored 
using a one-dimensional bulk model of the mixed layer originally pro- 
posed by Garwood (1976,1977) and expanded here to include a prescribed. 
vertical velocity. 

The first hypothesis to be examined is can alteration of the 
existing mixed layer model be reasonably done? Can a relatively 
simple mathematical transformation be applied to include a prescribed 
vertical motion without altering the basic physics of the entrainment 
process? The study will show that this in fact can be done, as revealed 
by a careful examination of the interactions at each step of the alter- 
ation of the original model. 

A second hypothesis to be examined is that vertical motion produced 
by long period internal gravity waves for which the Rossby number is 
small, through interaction with cyclical boundary conditions (wind 
stress and heat flux) can produce two-dimensional (x-z) structures in 
the temperature field of a simple infinite ocean having horizontally 
homogeneous surface boundary conditions. An objective will be to 
demonstrate the frequency dependence of the internal vertical motion 
upon the phases of the cyclical heat fields, and the wind stress. It 
will be shown that there is a non-linear interaction between super- 
imposed diurnal and semi-diurnal internal tidal waves. Finally, this 
Study will show the relative strengths of the resultant two-dimensional 
structures, and relate these to their importance in forming an improved 


thermodynamic ocean model. 





IT. MIXED LAYER MODEL AND VERTICAL MOTION THEORY 


A. MIXED LAYER MODEL 

To study the effects of internal vertical motion on mixed layer 
evolution, the Garwood (1976, 1977) mixed layer model was used in a one- 
dimensional form. Garwood developed an ocean mixed layer model using 
the Navier-Stokes equation of motion with the geostrophic component 
eliminated, the continuity equation for incompressible water, the heat 
equation from the first law of thermodynamics, the conservation of salt 
equation, and a linearized equation of state. The entrainment hypothe- 
Sis depends upon the relative distribution of turbulent energy between 
horizontal and vertical components and is offered as a plausible 
mechanism for governing both entrainment and layer retreat. 

Garwood suggested that planetary rotation influences the dissipa- 
tion of turbulence for deeper mixed layers and enables a cyclical steady 
State to exist on an annual basis. Furthermore, the rate of entrainment 
for the stable regime cannot be a simple linear extrapolation of the 
unstable situations. Unlike the atmospheric boundary case, most of the 
solar radiation does not penetrate the layer. Therefore, downward 
turbulent heat flux in the oceanic boundary layer is as important as 
the upward flux during the course of both diurnal and annual cycles. 
The non-linearity of the interface entrainment tendency parameter, 
which is greatest for stable surface boundary conditions, results in 
a modulation of the long term trend of the mixed layer depth by the 
diurnal component of surface heat flux. In this model, buoyant pro- 
duction is somewhat more efficient than shear production as a source 
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of energy for vertical mixing because of its unique effect on the 
vertical component of the turbulent velocity. 

The governing heat equation for this model assumes two forms: 
one for the mixed layer, the other for below the mixed layer. For 


the mixed layer, the equation has the form 


— —— 0 
Qg<T> = (wT), - (wT) + agi] Qdz 
at = pcyh -h (1) 


where < > denotes the mixed layer average and Q is the net heat flux, 
a product of downward solar radiation, upward turbulent heat flux to 
the surface, and surface back radiation. Below the mixed layer the 


heat equation becomes 








aT 
Ot Oz pC, az° (2) 


where K is the eddy viscosity, and -&wT)/dz is the term added for the 
vertical motion. 


The buoyancy flux at the entrainment zone is modelled according to 


ag(wT), = —<ww>'/2<E> 3) 
h 


which is derived from the local turbulent kinetic energy (TKE) budget 
where 

z= <-h (3a) 
In the entrainment zone, the buoyant damping attributed to (wt), is 
assumed to be balanced by the convergence of TKE flux. This convergence 
is modeled as a function of the TKE kinetic energy components and the 
vertical distance over which this energy must be transported. <w'w'> 
and <E> are the bulk values of the vertical component and the total TKE, 
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respectively. The prognostic equations for these mixed=-layer turbu- 


lence variables are derived using the bulk second-order closure methods 


of Garwood (1977), and are 


a (n<é>) = muy - agniwT), + agh[(wT), + WT), | 
a1’ 2 2k came 
- (<E>!2 + th)<E> (4) 

@ (h<ww>) = agh| (wT')_ > (w'T') + 
2 i 3 

(<E>-3<ww>)<e>"/2 - (<Esl@4fhyce> (5) 

3 

where 


Ts (r/p, ) 1/2 
with t being the magnitude of the surface wind stress and Py is the 
air density, assumed to be a constant, and 

R® = aghAT 

(Au? + Av?} 
is the bulk Richerdson number. The change in the total TKE, “E?, is a 
function of wind stress shear production, entrainment shear production, 
bouyant production/damping, and dissipation. The change in the vertical 
component, <w'w's, is a function of bouyant production/damping, pressure 
redistribution, and dissipation. The quantities AT, Au, and Av are the 
jump values between the mixed layer (assumed to be vertically uniform), 
and the underlying stable waters. These bulk TKE equations are solved 


algebraically assuming quasi-steady state for the TKE budgets 
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8 (h<E>) = 6) 
Ot 


and 


@(h<w w'>) = 0 
Ot 


With this assumption the computed values of <w'w'> and <E> are in 
phase with the surface boundary conditions--the wind stress, pus, and 
the buoyancy flux, ag(w'T’)), 

Once the entrainment bouyancy flux is known from (3), the downward 


fiuxes of momentum associated with entrainment at the base of the mixed 


layer, (wTuT)_) and (wivT)_. are given by 


-wu), = wAu = -(wTLA 
We eh hae ~ “hoe (6) 
-(wv) = wdAv = -(wT), Av 


Here the entrainment velocity is 


Ww = -(wT) 
; im (3) 


and the entrainment velocity is also defined as 





Ww = dh + w_ 
e ' h (9) 


if w, iS greater than zero. If Wo is not greater than zero, equations 
(10) and (11) are applied. 
1. 3oundary Layer Deepening by Entrainment 
The entrainment heat flux (W'TT) is determined algebraically 
from (3), and the steady-state forms of (4) and (5) in terms of given 
values of uy, (wT), h and Rj, This heat flux is then applied to 


the given temperature profile. The new h and new momentum and temperature 
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profiles are unique solutions provided that (i) the mixed layer is 
homogeneous, (ii) the value of h is just large enough to prevent 
unstable density profiles (@1T/@z greater than zero), and (iii) heat 
is conserved, 
2, Boundary Layer Shallowing 

Shallowing occurs when there is adequate vertical turbulent 
energy to transport heat down to the base of the existing mixed layer. 
The depth at which the vertical heat flux (w'T") vanishes establishes 
the new layer boundary, and will be the base of the new mixed layer. 
In this model, this occurs when (w'w') approaches zero. Then (3) is 


no longer applicable, and the steady state forms of (4) and (5) 


reduce to 
3 2 — TT, 
— mu? — (<E>@+fh)<E> + agh(wT ) 
’ “w En (10) 
6 =« <E>™ - (<E>'2 + fh) <E> + agh(w T )o 
2 (11) 


Without the planetary rotation term, or if the downward surface buoyancy 


flux is dominant, h is proportional to the Obukhov length scale 


tC - © ag(wT'), 


If, on the other hand, ow), = 0, h is proportional to u,/f. In 
general, however, the depth of the shallowing mixed layer is a function 
of the two non-dimensional parameters, h/L and hf/u,. 

Heat and momentum are conserved in the profile specifications 
and, because the mechanical energy budget requires conservation of 
potential energy on the time scale of the turbulence, a numerical 


procedure is used that also preserves potential energy. This 15S 
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intended to ensure that the deepening rate is as correct as possible 


when the mixed layer deepens to its former position. 


B. VERTICAL MOTION THEORY 

By following a water particle as it is vertically advected by a 
9rescribed velocity field, w(x,z,t), the distortion of any temperature 
profile, T(z), after a small, but finite length of time, (At), may be 
computed by non-linear theory. 


If we assume 


w = —Azcos(wt + kx) 


ay (12) 
is prescribed, then 
dz = —Azcos(wt + kx) 
dt D | (13) 


This may be reformed and integrated 


z to 
= -~A cos(wt + kx) dt 
[i = |- 


where a water particle starting at time t; at position z, is advected 

to position Zo at time to. Or following any constant isotherm 
initially at position 2], to Zo at ty gives 

i = z, exp{_A| sin(wt, + kx) — sin(wt, + kx)| } - 

wD (14) 


_ : =| 
ir t= ty t}<<w » then 


22 = exp {aat[cos(u(t +i) + wx] } (15) 


or 


z, 8 (16) 


ny 
UR 
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which is the general form approximated by the progression of equations 
in this study of internal vertical motion. Internal vertical motion 
is applied to the mixed layer model through equation (10). dh/dt, 

due to the internal vertical advection, Win? is derived from equation 


(16). 


For future study, the salinity budget can be handled in a like 


manner. 


23 





IIIT. EXPERIMENTS 


A, PRODUCTS 
The products used to identify features in the model runs and used 
for comparison are: 
1, Table (Sample Figure 1a) showing the hourly change in the mixed 
layer primary variables: 
H is the mixed layer depth in meters. 
T is the mixed layer temperature in degrees centigrade. 
S is the mixed layer salinity in parts per thousand. 
TAU is the surface wind stress in dynes per square centimeter. 


QNET is the net heat flux out of the mixed layer in watts 
per square meter. 


WOD is the value of the vertical motion parameter affecting 
the model in centimeters per second. 


DH is the change in mixed layer depth due to the vertical 
motion for that time period. 


2. Contours of the temperature variation, T(z,t), (Sample Figure 1) 
between z © 0 and z * -60 meters over the initial 24 hours of the model 
run. Varying resolution is used to better resolve the evaluation of 
T(z,t). Temperature contours are in degrees centigrade. 

3. Contours of temperature variation similar to 2, in a quasi two- 
dimensional mode. This is accomplished by running the model 24 times, 
each time passing the prescribed vertical motion through the 24 hour 
cycle of the model at a different, progressive phase. The 24 runs are 
then sorted by hour giving 24 data sets of 24 independent points with 
a spatial separation of kx from hour one through hour twenty-four. 
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DAY AR I 1g 3 TAU wNET 4DD Du 


14.956 34.9000 O50 4 o> 
14.339 34.900 0.$0 -100.8 
14.335 434.000 0.50 -184.3 
14.939 34.900 0.50 -248.4 
14.349 34.000 0.50 -288.7 
14.3958 34.000 0.50 <-302.5 
14.969 34.000 0.50 -288.9 
14.976 34.000 0.50 -248.7 
14.978 34.000 0.50 ~184.7 
100—os- 10.4 14.974 34.000 0.50 -102.3 
11 #10.7 14.963 434.000 0.50 -+ 4.1 
1.1 14.943 34.000 0.50 100.2 
13) aL 4 14.926 34.000 0.350 100.8 
14 11.8 14.909 34.000 0.50 100.8 
LE MIR 25 ci U 14.893 34.000 2.50 100.8 
2 12.4 14.377 34.000 0.50 100.8 
A za 14.862 34.000 0.50 100.8 
18 13.9 14.847 34.000 0.50 100.8 
190 13.3 14.833 34.000 2.50 100.8 
20 13.6 14.319 34.900 0.50 100.8 
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22 13.9 14.806 34.000 0.50 100.8 o. J as 13.882 

Pemelasto 14.793 34,000 9.50  #100.a 0. : a 

23. 14.4 14.781 34.000 0.50 100.8 oOo. : = 

Sails. ?)) 14.767 34.000 9.50 100.8 o. ; o. 

pee teeteso) 14,762 34.000 0.50 - 3.8 96.0 9.0 = 
2 2 18.1 214.762 34.000 9.50 -100.8 0.0 0.0 

2 3 18.2 14.768 34.000 0.50 -184.3 2.0 0.0 13.520 

> 4 5.3 14.777 34.000 0.$0 -248.4 0.0 9.0 
pecur is)) «14. ?@@) 34.000 0.50 -288.7 0.0 0.0 
2 6 18.8 14.901 34.000 0.50 -302.5 0.0 0.0 
2 7 15.4 14.914 34.000 0.50 -288.9 0.0 0.0 
Boaueises iaee23 34.000 0.50 -248.7 0.0 0.0 
2 9 15.6 16.829 34.000 0.50 -1846.7 0.0 0.0 
2 10 18.7 16.830 34.000 0.50 -101.3 9.0 0.0 
Sot? 5.3 14.825 34.000 0.50 - 4.1 6.0 0.0 
2 12 16.1 14.814 34.000 0.50 190.2 0.0 0.0 
2 13 16.3 14.803 34.000 0.50 100.8 2.0 0.0 
2 14 16.5 16.793 36.000 0.50 100.8 0.0 0.0 
2 1S 16.6 16.783 34.000 0.50 100.8 0.0 9.0 
>) femeieca 14.773 34.000 0.50 100.8 0.0 0.0 
Sue to” «C4. 763 «234.000 6.50 «100.8 0.0 0.0 
2 18 17.2 £4.1146.753 34.000 0.50 100.8 0.0 0.0 
Pete Slee) 14.743 34.900 0.$0 100.8 90.0 0.0 
2 20 17.5 16.734 34.000 0.50 100.8 9.0 0.0 
2 21 417.7 «214.725 34.900 0.80 100.8 0.0 0.0 
2 22 4217.9 16.726 34.000 0.50 100.8 93.0 0.0 
2 23 18.1 16.707 34.000 0.50 100.8 9.0 9.0 
2 26 18.2 16.697 34.000 0.50 100.8 0.0 0.0 


(a) 
= (b) 


y,) SE7E8) TEMP. LATER DEPTH 


14.9560 6.2613 t 
24.9360 7.2640 
14.9340 9.0113 
2 | 14.9399 9.4966 
z 4.3460 3.9230 
- 24.9580 3.2851 
= 14.9840 3.5392 (Ee ma | 
14.9780 3.8356 7a | 
14.9780 10.1260 ee soll 
24.9740 10.4177 ee 
S 14.9630 10. 7468 _ 
=] 4.9430 11.1046 
=| 34. 9280 11.9384 
2) :4.3089 11.7720 
-+ 24.8929 12.1007 
14.8770 12.4028 
14.8620 “2.7107 
14. 3460 13.9222 
=) 14,9320 13.2981 
=) 14.4490 13.5796 
q 34.8950 13.8668 
Sd ie. 7330 14.1371 
; 14.7400 1.3931 
14. 7670 14.6536 


Figure Ll. {a) Table of hourly values of 
problem parameters. H is the entrainment 

depth in m. T is mixed layer temperature 

in °C. S is mixed layer salinity in o/oo., 

TAU is the surface wind stress in dygjes/cm’. 
QNET is the net heat flux in watts/m*. WDD 

is the vertical motion parameter. DH is the 
hourly change in mixed layer depth due to 

the vertical motion in cm. (b) The contours 
of T9z,t) as evolved over a 24 hour period to 

a depth of 60 meters. (¢c) Successive temperature 
profiles from hour one through hour twenty four. 
The first profile is true to the temperature 
scale on the graph. Succeeding profiles are 
offset. 
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It is assumed that there is a two-dimensional internal wave in x and 
Z passing through this ocean system. Heat flux and wind stress are 
taken to be the same at all points, with no interaction between any 
of the points (advection is limited to the vertical). 

4. Progressive temperature profiles, T(z), to a depth of 60 meters 
(Sample Figure Ic). Only the first profile matches the horizontal 
temperature scale which is degrees centigrade. The remaining profiles 
have been consecutively offset. The accompanying table has the mixed 
layer temperatures in degrees centigrade, and the mixed layer depth 
in meters, 

9. Progressive temperature profiles as in 4, but in the two- 
dimensional mode of 3. Profiles are offset from the first, again, to 


provide an insight to the two-dimensional nature to the data. 


B. STANDARDS 

Prior to the examination of the addition of vertical motion to the 
mixed layer model (MLM), it is beneficial, here, to examine the results 
of running the mixed layer model in the two modes which will be used 
as the standards for comparison with the vertical motion modelling 
experiments. Both standard cases are runs of the mixed layer model 
from which the effects of evaporation and precipitation have been 
removed. This has been done to aid in the interpretation of later 
results by limiting the number of interaction variables to a manageable 
quantity. The first standard is a run of the MLM in the primarily wind- 
driven deepening mode. The second standard is for the MLM in a cycli- 


cal shallowing/deepening response to a strong downward heat flux. 
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1. Standard 1 (S1) 
This run was initialized with a temperature profile having a 
surface and mixed layer temperature of 15.0°C. At 5.0m the initial 


temperature decreases linearly with the form 


T(z) = T, - 5.0 x 1074*(Z-H) 417) 


Where Z 1s the absolute value of z in centimeters, and H is the mixed 
layer depth in centimeters. This profile extends to 200m. The result- 
ing profile is shown in Figure 2. 

The surface wind stress (TAU) is a constant 0.50 dynes per 
Square centimeter for the entire run. Net heating reaches a maximum 
of 302.5 watts per square meter at hour six (local noon). Net heat 
loss after sunset is a constant 100.8 watts per square meter. 

With no vertical motion the layer deepens continuously over 
the two-day period of the model run (Figure 1), as expected for pri- 
marily wind-driven mixing. In examining the first twenty-four hour 
period it is seen that the mixed layer deepens from an initial depth 
of 4.5 m to a depth of 14.7 m. The layer cools from 15°C to 14.767°C 
over the period, with a short period of mixed layer temperature 
increase of six hours related to that time when the net heat flux 
is sufficient to overcome the cooling attributed to the entrainment 
of cooler water into the layer. 

The salinity in the layer remains constant because there is 
no precipitation or evaporation allowed which could change the surface 
salinity, and the below layer salinity is assumed to be a constant 
34.00 part per thousand. Therefore, there are no sources of water 


with a salinity other than 34.00 parts per thousand. 
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Figure 2, Initial temperature and salinity orofiles for Case Sl. 
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The downward solar radiation, Qo increases from sunrise (hour 1) 
to noon (hour 6), and then decreases to zero at sunset (hour 12.) A 
constant is assumed for the net upward turbulent heat flux to the atmos- 
phere, plus back radiation of 100.8 watts per square meter for the 
entire day. 

The one-day response to this cycle can be seen in the contours 
(Figure 1b), The mixed layer deepens continuously over the first 24 
nours. B8elow the layer temperature, at any given depth not penetrated 
as the layer deepens, remains constant over time. In the layer the 
temperature is seen to cool slightly initially, then warm to a maximum 
value at about hour nine (1500 local time), and then cool again for 
the remainder of the period. 

In the T(z) graph this Beedtince is seen again (Figure Ic). 
Here it is seen that all of the slopes of the profiles below the layer 
are the same from hour one to twenty-four, indicating that with no 
vertical motion, the characteristics below the layer will not change 
with time, The included table shows the relationship of the temper- 
ature, in degrees centigrade, and the layer depth, in meters, and 
indicates the rate and magnitude of the deepening. 

The rate of deepening (Table 1) reflects the surface boundary 
conditions at any given time as well as layer depth and stratification 
below. As the run begins, the layer deepens quickly in response to 
the surface wind stress and the short distance over which it must act, 
as seen in equations (3) and (4). At hour one the layer deepens 
beginning at a shallow value where the wind flux has a high impact 


from equation (4), during a period where the heat flux into the layer 


vas, 








is a minimum, As the heat flux increases and the layer deepens the 

wind stress is relatively less effective in deepening, and the entrain- 
ment now effected by equation (5) slows until hour nine when the solar 
radiation goes to zero, and the deepening rate increases. This increase 
is relatively small because of the relatively large temperature jump 
(AT) in the entrainment zone by this time. After hour eleven (1800 
local time) the minor heat loss, small wind stress and increased depth 


involved all act to slowly decrease the deepening rate. 


Table 1. Deepening Rate in Meters Per Hour for Case S!] 


Hour Deepening Hour Deepening 
Rate Rate 
1 i613 ie 0.3296 
Z 1.0227 14 GF3370 
S 0.7271 15 0.3287 
a 0.4877 16 Gesu)? 
S 0.4242 17 0.3083 
6 O322) 18 OnSiiiS 
eo 0.294] 19 OR TASS 
8 0.2964 20 Oe2G lo 
9 0.2904 21 0.2872 
10 O.2on7 22 022703 
1] Ge3291 23 0.2560 
12 0.3580 24 0.2607 
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2. Standard 2 (S2) 

This run is begun with the same initial layer temperature as 
for standard 1. Here, however, the initial mixed layer depth is taken 
to be 49.5 meters. At 50.0 meters, the profile again assumes the 
linear variability of equation (18), extending to 200.0 meters 
(Figure 3). 

All other parameters are the same as those of S1 except for 


the prescribed surface heat flux. The net neat flux has increased to 
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Figure 3. Initial temperature and salinity profiles for Case Sd. 
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a maximum of 504.2 watts per square meter by assuming an extremely 
strong solar heating maximum (Q.) of 600 watts per square meter at 
local noon, 

In this case the layer depth remains at 49.5 meters for the 
first hour (Figure 21a) because the wind stress is not strong enough, 
coupled with the available heat energy, to cause further entrainment 
at this depth, and the solar heating is not yet strong enough to cause 
(w'w’) to go to zero so that saudtiens (10) and (11) are applied. By 
the second hour the net downward surface heat flux and the surface 
wind stress are sufficient to cause (w'w’) to go to zero, and a new h 
is formulated through the solution of equations (10) and (11), forming 
a new surface mixed layer with a depth of about 18 m. The layer then 
continues to shallow until hour 6 (local noon) in response to the 
increased downward heat flux. Stirring by the wind stress 1s relatively 
small as its effect is overcome by the buoyancy flux associated with 
net surface heat flux. The layer temperature increases as long as 
there is heating sufficient to overcome the minor cooling due to the 
entrainment of cooler below layer water which resumes at local noon, 
but is not strong until hour 10 (1600 local time). At that time the 
layer cools due to the combined effects of the surface heat loss, and 
the entrainment of the cooler below layer waters. 

The salinity, again, remains constant because there is no 
source of water of a salinity other than 34.00 parts per thousand. 

As shown by the contours (Figure 11b) the original layer of 50 meters 
shallows to hour 6, not forming a distinct layer, since equations 
(10) and (11) are applied at each hourly step, until hour 7. The 
transient thermocline prior to hour 7 is weak. This is seen more 


32 





clearly, perhaps, on the T(z) profiles (Figure llc). The originally 
sharp interface at 50 meters becomes diffuse with time due to small 
(Ky = 0.2 cm/s) below-layer diffusion. This is clearly seen in the 
T(z) profiles, and it is the reason for the non-horizontal contours 
near the 50 m depth. The variability of the temperature structure 
below the newly formed surface mixed layer is a function of the shal- 
lowing adjustment by equations (10) and (11) as applied from hour two 
through six. 

Again it is evident in the T(z) profiles that the profiles 
keep the same slope below the layer except for the originally sharp 
interface at 50 meters where diffusion of the layer occurs with time. 
This diffusion area is also seen to hold at a relatively constant 
depth with time, 

The T(z,t) evolution reflects the turbulent kinetic energy 
budget for this problem (Table 2). The initial decrease in h from a 
deep layer to a shallow layer is clearly the result of the surface 
heating and the surface wind stress acting to form a new, warmer, 
Shallower surface layer. Energy is increased in the surface mixed 
layer due to excessive solar insolation during the day. The Obukhov 
length is reduced from hour two through six, and the newly forming 
mixed layer adjusts accordingly. By hour seven the solar insolation 
begins to reduce. As this happens, the Obukhov length increases, 
allowing the newly formed layer to deepen. Deepening does not keep 
up with the Obukhov length because of the damping due to the energy 
transport which must take place during the mixing process. By hour 


11 the net heat flux is negligible decreasing to a constant surface 
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loss after sunset. The relative magnitude of the depth increase 
(@h/dt) with time when compared with that of $1 is due to the strength 
of the temperature jump (AT) at the layer, and the available energy 

in the layer. In Sl] a AT of the order of .4°C was encountered, 

while in S2, due to the slow formation of the new surface layer, a 

AT of the order .08°C is the maximum encountered. This allows for 
more rapid deepening. The rate of deepening increases after nour 16 
(2200 local time) because, as seen in the T(z) profiles (Figure 11c) 
the AT at the entrainment zone is decreasing with time due to the net 


upward surface heat flux. 


Table 2. Deepening Rate in Meters Per Hour for Case S2. 


Hour Deepening Hour Deepening 
Rate Rate 

] 0.0152 13 0.8464 
G -31.4079 14 0.3317 
3 -1,9247 15 0.8252 
4 -3.1646 16 0.8235 
5 -0.9575 17 0.8278 
6 -0.2935 18 0.8404 
7 0.3766 19 0.8596 
8 0.3445 20 0.8879 
9 0.4944 2] 0.9186 
10 Os5see 22 0.9516 
17] 0.6793 £3 0.93870 
2 0.8371 24 1.0256 
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C. EXNPERTMENTS 

The addition here to the existing theory of Garwood (1976-77) is 
approached in a step-by-step procedure, gradually increasing the com- 
plexity of the problem. <A quick reference for this progression can 


be seen in Table 3. A basic assumption that was made in all of the 
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Table 3. 


Initial conditions, 


boundary conditions, and progression 


of vertical motion inputs for al] experiments. 


For all Cases: 


Case 


Sl 
$2 
la 
Ib 
Ic 
Id 
Ila 
IIb 
IIc 
IId 
Illa 
IIIb 


IVa 


IVb 


Va 


Net surface heat loss is constant at 100.8 watts/m?. 


Surface wind stress is constant at 0.5 dynes/cm2. 


The initial layer temperature is 15.0°C. 


The salinity is constant to 200 m at 34.00°/oo. 


No precipitation or evaporation is allowed during the run. 


The vertical motion is exponentially decreasing to the surface. Maximum motion 
is held at 36.0 cm/hr at 10 m, except for Case IVb, where it is 18.0 cm/hr at 


10 a. 


B. Cc. 
Q, max. 
(watts/m) 
400.00 


600.00 


Sl 
S2 
$1 
$2 
Sl 
S2 
S1 
S2 
$1 
S2 


S2 


S2 


$2 


S2 


iC. 
Temp. MLD 
Profile (m) 
Fig. 2 4.5 
Fig. 3 49.5 
S1 S1 
$2 S2 
SI S] 
S2 S2 
$1 $1 
$2 $2 
$1 Sl 
S2 $2 
S1 S1 
$2 oy. 
$2 S2 
S2 S2 
S2 Sz 
Sz S2 


Vertical Motion 


None 


"None 


Constant downwelling during model run. 


Constant upwelling during model run. 


Pulsed downwelling, 12 hour cycle. Phase 
shifted to peak at hours 5S, 18, and 24. 


Pulsed upwelling, 12 hour cycle. Phase 
shifted to peak at hours 5S, 18, and 24. 


Single diurnal wave passage. Phase shifted 
to give maximum downward motion at all hours. 
Hours 1, 7, 13, and 19 discussed. 


Continuous diurnal wave system passage. Same 
phase shifts as Case III. Hours 6, 12, 18, 
and 24 discussed. 


Continuous semi-diurnal wave system passage. 

Phase shifted as in Case III. Runs with peak 
motions at hours 6 and 18, and 12 and 24 are 

discussed. 


Superimposed diurnal and semi-diurnal continuous 
wave system passage. Both waves in phase. 
Input phase shifted as in Case III. 


As in Case Va, but phases of diurnal and semi- 
diurnal waves varied relative to one another. 
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work was that the vertical motion moved the water up and down, but did 
not directly induce turbulence (i.e., mixing). This is a valid 
assumption if no shear is associated with this motion. 

The cases of progressively complicated vertical motion examined 
were: 


1, Linearly (in z) varying vertical motion, constant over the 
entire period. 


a. Upwelling 
b. Downwelling 

2. Linearly varying vertical motion pulsed in time. 
a. Short daily upwelling/downwelling events. 

3. Linear internal wave motion. 


a, Passage of a single wave pulse through the experimental 
point over the 24 hour daily cycle. 


4, Non-linear wave variations coupled with diurnal and semi- 
diurnal boundary conditions (heating, wind). 


5. Non-linear response to’ superimposed waves of different periods. 
Diurnal and semi-diurnal tidal waves are combined and phase- 
shifted to determine non-linear interactions. 

le Gage I 
Both upwelling and downwelling are realistically induced with 

the realization that in a mixed homogeneous layer the vertical dis- 


placement due to a linear w(z) can be shown to be exponential, going 


to zero at the surface. Mathematically we know that 


ah = w(z) (18) 
Wi 
w(z) = AZ 
D 


where D is a Scaling depth, it follows that 
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Qa 
= 
or 


or 


dz = Ad 
z D 
Integration from t, to ty yields 


2. = z; eat 


which is the equation upon which the following work is based, and which 


(19) 


agrees with equation (16) derived earlier. 

The temperature and salinity profiles must reflect “correct" 
variation if equations (16) and (19) are correct. Upwelling and/or 
downwelling will not cause isolines to break the surface or each other. 
Mixing must be done to induce this. Therefore, above the mixed layer 


depth Te and ae must equal eae and Sold when there is vertical 


wW 
motion without mixing or surface exchange. Below the entrainment depth 
the profiles must be altered to reflect the type of vertical motion 
(upward or downward) in accordance with equation (16). Figure 4 shows 
the expected result from upwelling on a temperature profile. The mixed 
layer depth nas decreased from time one to time two (after upwelling), 
and the slope of the profile below the mixed layer has been decreased. 
In all the cases the AT at the entrainment zone has been preserved. 
a. Constant Downwelling (Case Ia) (Figure 12) 

The initial conditions used are the same as those for S|] 
(See Table 3). A constant downward vertical motion is introduced giving 
the vertical motion parameter, WDD, (WDD = A t/D) a value of .03600, 
which equates to a value of vertical motion of .01] cm/s or 36 cm/hr at 
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ten meters. This is a rather large value, but it demonstrates the inter- 
actions well. The vertical motion is applied beginning the first hour, 
and 1s not applied as the initial condition on the temperature and 
salinity profiles. 

The result of this downward motion is, as expected, a more 
rapid deepening of the surface mixed layer (Figure 12a) when compared 
to SI. Comparing the depth change due to the mixing in this case (Table 
4) to the S1l mixing, the same sequence of events as described for $1 is 
seen, and the values are comparable, indicating that the physics of the 
model has not been altered by the addition of the vertical motion para- 
meter in the downwelling mode. It is noted that the effect of the 
vertical motion quickly becomes greater than the effect of mixing due 
to the large value of vertical mixing assigned in this run of the model. 


Table 4, Deepening Rate in Meters Per Hour for the 
Mixing and Vertical Motion of Case Ia. 


Mixing Vertical Motion Mixing Verticle lotion 
Hour Deepening Deepening Hour Deepening Deepening 
Rate Rate Rate Rate 
1 L.7613 0.2295 13 0.3201 Qeo72o 
2 0.9705 O0:2735 14 03169 0.6055 
5 0.6006 Ue3077 iis 0.3134 0.6392 
4 0.4759 0.3365 16 0.3102 0.6740 
5 Oweabie5 0.3613 17 0.3068 7099 
6 0.2945 0.3858 18 0.3034 0.7471) 
/ 0.2444 0.4089 19 0.3017 0.4850 
3 0.2282 0.4323 20 0.2963 0.8252 
9 0.2394 0.4569 2] 0.2949 0.8663 
10 0.2513 0.4828 22 0.2914 0.9087 
7] 0.2793 0.5108 23 0.288] 0.9526 
12 003219 0.5413 24 0, 2094 eogs2 
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The sequence is vividly displayed in the contours (Figure 12b) 


where it is seen that the temperature contours deepen as time progresses 
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as one would expect in a downwelling case. The surface heating due 

to the incoming solar radiation is noted again, followed by cooling 

as the sun sets and heat is then lost at the surface and colder water 
is entrained at the bottom of the mixed layer. It is also noted that 
in only twenty-four hours, the vertical motion added here has deepened 
the layer to almost twice the value of S1. 

The T(z) plot (Figure 12c) demonstrates the deepening 
involved and the effects on the temperature profile below the mixed 
layer. Note the increased slope (dT/dz) with time, as warmer near- 
surface temperatures (water) are carried downward by the downwellins. 
Also, one can see the effect of the preservation of the AT ipnetiie at 
the entrainment zone from the mixing step through the vertical motion 
progression in the model where at times there is an apparent discon- 
tinuity in the profile. This does not appear to cause a problem in 
the execution of the model, and will be seen in many of the following 
T(z) plots. One can further see that the vertical motion addition has 
caused a slight decrease in the jump strength of AT at the layer over 
that in Sl. This can be attributed to the reduced mixing for the 
larger mixed layer depths. Where, with no downwelling a AT of 1.0°C 
may have occurred in 18 meters, it now occurs at a greater depth. 
This effect is transmitted into the entrainment zone as this below 
layer water is mixed in, causing the AT to be reduced. 

b. Constant Downwelling (Case Ib) (Figure 13) 

The initial conditions are those of S2 (See Table 3). A 

constant downward motion is introduced equal to that used in the 


Drevious case. 
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The result of the downward motion is again, as expected, 
a more rapid deepening of the layer after the new surface mixed layer 
is formed through the surface wind stress and heat flux as applied by 
equations (10) and (11) (Figure 13a). Comparing the depth change due 
to the mixing here (Table 5) to the S2 mixing shows that the values 
are of the same order, which shows that the energy relations of the 
original model are not impaired by the addition of vertical motion 
in this mode of the model either. Here it is noted that mixing drives 
the production of a new surface mixed layer beginning at hour 2 (0800 
local) (Figure 13b). The shallowing due to the surface heating is 
dominant until hour seven (1300 local) when the downwelling becomes 
the dominant factor until hour nine. From hour nine to hour fifteen 
the effects of mixing and downwelling are of the same order. After hour 


Table 5, Deepening rate in meters per nour for the 
Mixing and Vertical Motion in Case Ib. 





41xI1ng ertical ‘lotion Mixing Vertical Motion 
Hour Qeepening Deepening Hour Deepening Deepening 
Rate Rate Rate Rate 
] 0.0152 1.8150 iS: 0.8093 0.6986 
2 -33.4052 0.6571 14 0.7343 0.7511 
S =2.5202 0.5885 15 0.7445 0.8059 
4 -3.7245 0.4735 16 05713] 0.8616 
5 -1.4432 0.4392 ie 0.7029 0.9190 
6 - ./286 0.4286 18 0.6768 0.9775 
Z 0.1233 0.4489 19 0.6716 1.0379 
8 0.3275 0.4773 20 0.6628 1.1002 
9g 0.4266 0.5104 2] 0.6420 1.164] 
10 0.5737 0.5504 22 0.6316. 1.2299 
1] 0.6446 0.5942 23 0.6224 1.2978 
12 0.8003 0.6986 24 0.6120 1.3678 


rrr nee fm OE nS a eS 


fifteen the downwelling again, slowly, becomes the dominant factor for 


layer deepening. 
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It is seen that the layer temperature, again, immediately 
begins to warm (Figure 13c), and in fact, while the dominant factor is 
the large amount of surface heating, the layer warms at exactly the 
same rate as in S2, After hour ten the layer again begins to cool, but 
because of the deepening of the mixed layer due to downwelling, the 
mixing of the surface wind stress is reduced. This, along with the 
fact that a deeper layer has available more mass of water to dilute 
the effect of the entrained water, results in less cooling over the 
remaining twelve hours. 

Because of the downwelling effect on diffusing the temper- 
ature gradient it is also seen (Figure 13b) that the entrainment zone 
(layer boundary) remains rather diffuse until hour eight, while in S2 
it was sharpened by hour six, This helps explain the isagnitude of the 
deepening effect in the period from hour nine to fifteen, thougii the 
depth of the entrainment zone is increasing faster than that in S2 and 
reaches a difference of five meters from that in S2 by hour fifteen. 
The AT across the bottom of the layer for the downwelling case is 
Slightly less than that in S2, allowing the mixing to be a more impor- 
tant factor for deepening in this case. The following decrease in the 
mixing is then due to the ever increasing depth through which the mix- 
ing must take place. The wind stress must act over distances that are 
twice as large as those in S2, and although the AT at the layer depth 
is not as great, the effect of the depth is quite noticeable. 

The original deep layer is seen to begin to diffuse, but 
is quickly carried below the 60 meter level of the graphs. As stated 


above, the new surface mixed layer begins to show a significant AT 
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at hour eight. Extremely close examination of the T(z) profiles for S2 

and the profiles for this run does show a smaller AT at the entrainment 

zone for this case during the period from hour seven to hour sixteen. 

After hour 16 there is no discernable difference in the AT in the profiles. 
c. Constant Upwelling (Case Ic) (Figure 14) 

The initial conditions are those of $1. A constant upward 
vertical motion is introduced beginning at hour one giving the vertical 
motion parameter, WDD, a value of .03600, or about 36 cm/hr at ten meters. 
Again, this 1s high, but it is done with the realization that later work 
will be with wave motion where this will not be an unrealistic value, 
and using it here provides a good benchmark for the effects produced by 
this motion. 


Table 5, Deepening Rate in Meters Per Hour for the 
Mixing and Vertical motion of Case Ic. 


Mixing Vertical Motion Hour Mixing Vertical Motion 


Hour Meepening Deepening Deepening Deepening 
Rate Rate Rate Rate 
1 1.7613 -.2214 13 0.4239 -.3026 
2 vv. 1081 -.2526 14 0.4528 -.3079 
3 0.6746 -.2675 15 0.4845 -.3142 
4 0.6103 -.2/96 16 0.4774 -.3199 
5 0.4954 -.2873 7 0.2909 -.3189 
6 0.3417 - .2892 18 0.2871 -.3178 
7 0.3184 -.2902 19 0.3046 -.3173 
8 0.3097 -.2909 20 0.3223 -.3175 
9 0.3133 -.2917 (a 0.3407 -.3183 
10 0.3254 -.2929 22 0.2596 -.3163 
til 0.3435 -.2947 23 0.2739 -.3148 
12 0.3938 -.2982 24 0.2891 -.3138 





The result of the addition of this upward motion is, as 
anticipated, a slowing of the deepening rate of the mixed layer (Figure 


14a) when compared to Sl. In fact, due to the size of the vertical 
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motion parameter, an effect is produced which almost holds the mixed 
layer at a constant depth. This is displayed more dramatically after 
two model days have been executed, and can be seen to reach a near 
equilibrium condition after four to five days. Comparing the deepening 
rates due to the effect of the mixing in the model and the vertical 
motion in this case (Table 5) to the $1 mixing it can be seen quite 
clearly that the effects described for the mixed layer model standard, 
S1, are still present. From equations (1) and (3) the slightly 
increased magnitude of the cooling of the mixed layer by entrainment 
and the subsequent heating are attributed to the reduced depth over 
which the wind stress and heat flux must be distributed. Again, it 
can be seen, that the model TKE balance, equations (3-5), is not 
violated by the addition of this upwelling motion. 

The layer temperature initially falls (Figure 14a), then 
rises until 1400 local time. It then falls again in response to the 
mixing and the surface heating cycle, as in SI. The heating of the 
layer is not as great here as in the downwelling case (Figure 13), 
because the entrainment is sufficiently enlarged due to the smaller h, 
and because cooler water becomes available for entrainment more quickly 
due to the upwelling. This iS more clearly seen in comparing the 
contoured fields, T(z,t), (Figure 14b) of this case with those of the 
downwelling case (Figure 13b). By the end of the twenty-four hour 
period this has resulted in a very strong temperature jump, AT, at 
the entrainment zone, as is seen best in the T(z) profiles (Figure 14c). 
It is seen that AT grows with time, 

The strengthening of the jump, and the lifting of the 
cooler water is demonstrated clearly by the contours, and it is 
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clearly evident the contours not only lift, but also merge in time, 
indicating the AT grows with time. This effect is seen in the T(z) 
profiles by the increased gradient of the profiles with time. 

Equilibrium is reached in this case as the jump gets too 
large for the mixing to have an effect greater than the effect of 
the vertical motion, Equation (8) shows that dh/dt vanishes as Wo 
approaches wp. 

d. Constant Upwelling (Case Id) (Figure 15) 

The initial conditions applied are identical to those of 
S2, A constant upward vertical motion is introduced at hour one of 
about 36 cm/hr at ten meters, as before, 

The net result of the vertical motion after 24 hours is a 
reduced layer depth in comparison to the case having no vertical motion 
(S2). The total reduction in depth of six meters, nowever, is not as 
great as the previous case (Ic). After two days it is seen that the 
mixed layer iS approaching an equilibrium depth of near eleven meters 
(Figure 15a). 

The deepening rates (Table 6), when compared to those of 
$2, demonstrate that the basic model is still preserved in this run. 
From hours two through five the model is in the shallowing mode, and 
the new layer depth is established by equations (10) and (11), while 
there is a massive influx of solar radiation to the sea surface. The 
increased deepening rate due to the mixing seen from hour twelve 
through twenty-three is attributed to the effect of equations (1), 
(4) and (5). The mixed layer is maintained at a shallow level by the 
upward vertical motion, therefore, the mixing takes place over a much 


smaller distance than in S2. 
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Table 6, Deepening Rate in Meters Per Hour for the 
Mixing and Vertical Motion of Case Id. 


1xing ertical ifotion Hour Mixing ertical Motion 
Hour Deepening Deepening Deepening Deepening 
Rate Rate Rate Rate 
] 0.0352 -1,7509 13 0.9465 -.4751 
2 -29.4675 -0,6470 14 0.9470 -.4917 
3 -1,3404 -0.5767 15 1.0159 -.5103 
4 -2.6147 -0.4639 16 1.0611 -.5298 
5 -0.5166 -0.4292 17 1.1449 -.5515 
6 Glen -0.4186 18 1.2301 -.9755 
7 0.3085 -0.4148 19 1.3043 -.6013 
8 0.4578 -0.4163 20 1.4087 - .6298 
9 0.5885 -,4224 21 los 7 - .6626 
10 0.6775 -0,4314 22 1.6638 - .6980 
1] 0.7113 “0.4413 23 1.3743 -./219 
12 0.9078 -0.4578 24 0.6909 -./208 


Further, there is a concentration of energy in this shallower 
layer. The great variance in the last hour is a direct result of the 
magnitude of the upwelling, which, by hour four, has set up a large 
gradient in the temperature (Figure 15b). By hour twenty-four a gradient 
has been set up with such a large slope that mixing is becoming very 
difficult. Continued upward motion of this magnitude, without shal low- 
ing as is seen on the second day, would quickly make it extremely 
difficult to entrain much water into the mixed layer. 

The progression in the change of the mixed layer temperature 
(Figure 15a) shows that as the model is shallowing, governed by equation 
(10) and (11), the same adjustment takes place as is seen in S2. The 
temperatures produced are exactly those seen in S2, This is another 
indication that the model is performing as before. Since w, is zero, 
the upwelling has no effect on the setup of the new layer. If the H 
values had been listed immediately after the mixing phase of the compu- 


tation, this could be seen more clearly. After the layer again begins © 
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to mix downward, however, the upwelling begins to have some effect on 
T and h due to the upward transport of the entrainment zone. Though 
the layer is one to two meters shallower than in S2, the temperature 
of the layer is the same as that in the deeper standard to hour ten. 
At hour ten the heat flux of equation (1) has been reduced sufficiently 
so that wind mixing again becomes the dominant deepening force. At 
this point the upwelling has moved water cooler than that in S2 at the 
same time to the entrainment zone. The change in entrainment rate due 
to the upwelling, the entrainment of cooler waters, and the net surface 
heat loss after hour eleven combine to produce a cooling of .147°C from 
the maximum warm temperature, as compared to .109°C for S2. For the 
mass of water involved, and the relatively short time involved, this 
is a Significant difference. 

The effect of the upward vertical motion is perhaps best 
Shown in the T¥z) diagram (Figure 15c), where it is seen that the model 
initially adjusts, through equations (10) and (11), to a new surface 
layer which slowly strengthens, in terms of the jump, AT, as it en- 
trains. As the process continues the profiles below the layer rise 
by upwelling and the initial deep interface is smoothed by diffusion 
until it disappears. By hour twenty-four a sufficiently stable pro- 
file has been set up to make deepening of the layer more difficult. 
This indeed is the case in the second day of the run (Figure 15a) as 
the layer reaches a temporary equilibrium, It is temporary, in that 
once the profile has been reached which makes deepening sufficiently 
difficult through equations (4) and (5), the dominant force becomes 


the upwelling, w_, of equation (9). 
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The contours (Figure 15b) further demonstrate the effect 
of the upward lifting of the cooler water, as the contours tend to 
converge with time. The initial interface is quickly lost, as the 
model smooths this into its uplifted profile. 

2. Case II 
The next step of this progressive series of tests is to explore 
the effects of a simple pulsed motion varying in time. This is accom- 
plished by varying the input motion in a sinusoidal pattern. By using 


the equation 


w(t) = AAsin(27t + kx) 
(20) 


as the particular form of equation (12) applied for this case, where 
AA is the velocity scale factor such that AA = Az/D is scaled for an 
appropriate motion at 100 meters, and limiting the value of Wo(t) to 
only values in one direction (+ or =), and zero, a pattern of motion 
can be produced such as that seen in Figure 5. A pulse of vertical 
motion of twelve hour duration is produced, followed by a twelve hour 
period of no motion. To determine the relationship of this motion 
with the cyclical boundary conditions one need only shift the phase 
of the pulse relation. This will cause the motion to peak at varying 
selected times during any twenty-four hour period run. 

a. Pulsed Downwelling (Case Ila) (Figure 16, 17 and 18) 

The initial conditions are the same as those of Sl. Three 

different phase relationships of the downwelling to the cyclical 
boundary conditions are presented. In the first (Figure 16) the 


maximum of the downwelling occurs 180° out of phase with the heating 
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cycle (that period of the model defined as night, or the cooling cycle). 
In the second (Figure 17), the downwelling occurs during the last seven 
hours of night, and during the first five hours of daylight. The max- 
imum downwelling occurs six hours before the heating maximum. Finally, 
in the third case (Figure 18), the downwelling occurs in phase with 
the maximum heating cycle. In all three cases, the downwelling event 
is a twelve hour pulse with a maximum value of 36 cm/hr at ten meters. 
(1) 180° Qut of Phase (Figure 16). For this case it is 
seen that the model evolves exactly as does S1 until hour thirteen 
when downwelling begins (Figure 16a). All the fields produced are 
exactly the same as in S1 until hour thirteen for which the evolution 
has been discussed, 

After hour thirteen the downwelling begins to trans- 
port the surface water downward, and causes a decrease in the jump, AT, 
strength of a few hundredths of a degree. This combined effect causes 
the deepening of the mixed layer to increase here over that of Sl by 
about four meters over the twenty-four hour period. Additionally, it 
decreases the change of temperature of the mixed layer because the layer 
is initially deeper, thereby, reducing mixing, and increasing the dilu- 
tion of the entrained waters, and by causing entrainment of water which 
is not as cool as that of Sl, over the same period. These effects are 
apparent in the contours (Figure 16b), and in the T(z) profile (Figure 
16c). The downwelling effect is seen most dramatically in the contours, 
where the downward pulse is evident, along with the described temper- 
ature variations. 

(2) Waximum Downwelling Before Heating (Figure 17). Here, 
the downwelling is applied during the last six hours of the mixing, 


90 








and during the first five hours of the heating cycle. The six hours 
of downwelling during the last six hours of night serve only to deepen 
the layer in the manner previously described. The first five hours of 
downwelling coupled with the heating and wind stress through equation 
(9), produces a noted change in the temperature fields (Figure 17c) 
from that of SI. This coupling leads to an increase of the layer 
temperature at the peak heating period of .008°C for this run over 
that of the previous (8T increase), and over that of the standard, S1. 
After twenty-four hours this difference has decreased to .005°C as 
compared to the previous run, and by two days the difference in temp- 
eratures is negligible, though the layer depth is slightly shallower 
after two days for this run. By only partially coupling the effects 
of downwelling and heating, a shallower mixed layer has been produced 
with somewhat warmer temperatures. 

The surface temperature evolution and the below layer 
evolution are best revealed by the contoured fields (Figure 17b). It 
1S apparent here that another phenomenon is present. As the downwell- 
ing commences the near layer contours do in fact diverge as expected. 
The deeper contours, however, initially narrow as the pulse forms, 
then widen to the below layer temperature spreading effect described 
previously. This convergence during the peak motion period is a direct 
result of the mathematical treatment of the below layer temperature 
profiles through equation (16). 

(3) Coupled Heating and Downwelling (Figure 18). dere 
the downwelling peaks one hour before the maximum surface heat flux. 
The combination of downwelling and wind stress, equation (9) through 
equations (16) and (8),are sufficiently strong to cause the effect of 
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the heat flux, equation (1), to be negligible, and the layer deepens 
to almost three meters more than that of SI by hour eleven. Also, 
the layer temperature has been increased significantly over S1 due 

to the downwelling effect on the layer which causes the entrainment 
of less water than previously during the period of maximum surface 
heating. By combining the downwelling with the heating, a maximum 
temperature change of .211°C for SI has been reduced to .162°C for 
this run over the first twenty-four hours. This has also produced a 
mixed layer which is .06°C warmer and three meters deeper than S] 
after only twenty-four hours. This reflects a substantial change in 
the energetics of the system, and may be a significant mechanism for 
the generation of horizontal temperature structure. After the pulse, 
the layer is mixed by the action of the wind stress in the maceenes 
of a smaller AT--reduced by the downwelling accompanied by a smal] 
downward surface heat flux. This explains why the deepening is not 
as great after the downwelling pulse as would be expected if one con- 
Sidered only the reduced jump strength. 

In combination these three phased runs demonstrate 
the effect of the interaction of the downwelling and the diurnal 
heating cycle. In all cases a relative heating of the mixed layer 
is seen because of the depth increase. The deeper layer reduces the 
effect of the mixing, and reduces the entrainment. As the downwelling 
event couples with the heating, the final mixed layer depth is reduced, 
over that of the uncoupled cases, and greater surface temperatures 
(layer temperatures) are achieved. The rudiments of a pattern in 
surface features which could be induced over a large ocean by simply 
varying the coupling of these two basic effects iS now seen. 
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b. Pulsed Downwelling (Case IIb) (Figures 19, 20 and 21) 

The initial conditions are those of $2. The magnitudes 
of the downwelling, and the phase relations are the same as the 
previous case. 

(1) 180° Out of Phase with Heating (Figure 19). Because 
the downwelling event occurs after the surface heating has stopped (at 
night), the model has gone througn the shallowing mode, equations (10) 
and (11), and produced a new surface mixed layer exactly as in S2 dur- 
ing the daylight period. The effect of the downwelling is seen at hour 
fourteen when mixed layer deepening is first observed. The cumulative 
effect of the downwelling is to deepen the layer by about five meters 
more than in S2 in the twenty-four hour period (Figure 19a). The layer 
temperature, due to the downwelling effect, remains .011°C warmer after 
the period. Destruction of the original layer at 50 m is caused by 
diffusion, and the model's mathematical smoothing in the production of 
new temperature contours due to the downwelling in accordance with 
equation (16). 

(2) Maximum Downwelling Before Heating (Figure 20). The 
shallowing of the model through equations (10) and (11) during the 
Surface heating period overrides the downwelling, so that by hour five 
this run is essentially the same as S2, with only a slight temperature 
difference. Comparing the contours (Figure 11b and Figure 20b) shows 
that this run is the same as S2 until the downwelling pulse begins at 
hour nineteen, except the early downwelling has deepened the layer 
and moved the cooler water away from the entrainment zone, relative 
to the temperatures in SZ. The layer then deepens due to wind stress 
mixing, and downwelling. The net result of this coupled action is a 


4 


5S: 





layer about four meters deeper than S2 after twenty-four hours, but 
with a temperature essentially equal to that of S2. 

(3) Coupled Heating and Downwelling (Figure 21). dere 
the shallowing mode is still the dominant factor, with the downwelling 
playing only a minor role in that it deepens the newly forming layer, 
and moves the relatively cooler water away from the entrainment zone. 
Because the model is in the shallowing mode, aS We approaches Zero, 
the new value of h is dependent on equations (10) and (11), as is the 
temperature of the new mixed layer through equation (1). Since the 
downwelling adjustment is done after the evaluation of the shallowing 
h, the downwelling does deepen the newly formed layer, but no temper- 
ature adjustment is made. This results in a layer that appears to be 
deeper than in S2, and yet maintains the same temperature. The net 
result of these interactions after twenty-four hours is a layer that 
is one meter deeper than in S2, and, .012°C warmer. Here the coupling 
effect was negligible. 

In combination these three runs of the model in the 
downwelling mode with S2 initial conditions showed the opposite pro- 
gression of those runs using the S] initial conditions. The maximum 
variations of temperature and mixed layer depth occurs in the case of 
180° out of phase deameai ling with the heating cycle, where a deeper, 
warmer mixed layer is formed during the twenty-four hour period. In 
the case of the S1 initial conditions the maximum deviations occurred 
during the coupled run. 

C. Pulsed Upwelling (Case IIc) (Figure 22, 23 and 24) 


The initial conditions are those of Sl. The phase 
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relationships used are the same as in the pulsed downwelling case, as 
are the magnitudes of the effects. 

(1) 180° Out of Phase with Heating (Figure 22). The 
effect here is negligible until hour thirteen, and this run is exactly 
the same as S] until then. At hour thirteen the upwelling begins to 
increase the magnitude of AT. Colder water is moved up toward the 
entrainment zone below the layer, and the layer itself is moved upward, 
as the wind stress acts to deepen the layer by mixing. At hour twenty- 
one it is seen that the combination of the upwelling and the magnitude 
of AT are great enough to actually shallow the layer for that hour. 

The net effect of the upwelling in this case is to 
create a mixed layer that is about 2.6 meters shallower than in S]1 
over the twenty-four hour period, and .03°C cooler (Figure 22a). Also, 
there 1S a maximum 8T in the mixed layer of .241°C, where mincing only 
produced a 81 of .211°C in SI. This has the potential of producing 
an observable surface temperature effect. 

The effect of contour convergence during the pulse 
mentioned in the downwelling case is seen here again (Figure 22b). 
The overall effect of the upwelling, to converge the contours from 
before to after upwelling, iS seen as expected. 

(2) Maximum Upwelling Before Heating (Figure 23). In 
this case the upwelling immediately begins to lift the entrainment 
zone against the effect of the wind stress, causing an increase in 
the effect of the mixing, and reducing the water available for 
dilution of the entrained, cold water. In equation (9), w_, is of 


opposite sign to that of Wo» reducing, or changing the sign of dh/dt. 
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This results in cooler layer temperatures over the same period as 
described for Sl], and shallower mixed layer depths. Numerically, 

the net result of this effect over the twenty-four hour period is 

to produce a layer that is 2.5 m shallower than that of S1, and .04°C 
cooler (Figure 23c). The maximum 8T of the layer over the period is 
240°C, or about the same as in the previous upwelling case. 

(3) Coupled Heating and Upwelling (Figure 24). The 
coupling of the upwelling and the heating cycle produce a mixed layer 
that is only 1.9 meters shallower than that of S] in the twenty-four 
hour period, but a layer that is much cooler. The reduction of the 
layer depth during the heating period causes the maximum temperature 
to be .013°C colder, and causes the layer temperature at the end of 
the twenty-four hour period to be .071°C colder than SI, a significant 
change. The maximum 8T of the layer for this period is .269°C. 

As in the downwelling case with the shallow initial 
layer conditions, the maximum surface variations in temperature is in 
the coupled upwelling case. 

d. Pulsed Upwelling (Case IId) (Figure 25 and 26) 

The initial conditions are those of S2. The phase rela- 
tionships used are the same as those in the previous three model runs. 
Magnitudes of the upwelling are the same as in.the previous run. 

(1) 180° Out of Phase with Heating (Figure 25). The 
shallowing proceeds as in S2 for the first six hours, and a new layer 
is formed (Figure 25b). At hour thirteen the upwelling begins and 
several effects are noted. First, the deepening immediately slows, 
relative to that in S2, as the upwelling lifts the layer, and 
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increases the AT at the entrainment zone. Second, the layer temperature 
decreases faster relative to that in $2, due to the shallowing of the 
layer, and entrainment of relatively cooler water into the layer. 
Finally, there is an effect not previously observed. The initial layer, 
which is being diffused in the previous runs is now reinforced during 
the upwelling (Figure 25c), and a two-layer system is formed over the 
first twenty-four hours. 
The net result of this case is a surface layer that is 

2.4 meters shallower and .022°C cooler than S2 over the period. This 
leads to a slightly increased surface signature over the non-upwelling 
case. Because of the reinforcement of the initial layer there is also 
a potential for a substantially recognizable T(z) signature for this case. 

(2) Maximum Upwelling Before Heating (Figure 26). The shal- 
lowing and upwelling combine to make this a relatively uninteresting 
case. Though the large upwelling has some initial effect, the fact 
that w, goes to zero forcing the model to be driven by equations (10) 
and (11) soon overcomes the upwelling, and by hour four this case is 
identical to S2 and the previous run with the exception that the up- 
welling has advected cooler water closer to the entrainment zone. 
The difference in water temperature and mixed layer depth between 
this case and S2 are negligible by hour eighteen. At hour nineteen 
upwelling begins again, increasing the relative effect of the wind 
stress mixing, causing increased entrainment of cool water. The net 
effect of the upwelling and the wind stress mixing over the depth 
involved is to produce a layer that is 2.2 meters shallower than in 


S2, but at only .004°C cooler after twenty-four hours. This would 
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not produce a noticeable effect at the surface. The two-layer type 
system seen in the previous case could, again, provide a signature 
to this sort of upwelling event. 
(3) Coupled Heating and Upwelling (Figure 26), The 
Shallowing prescribed by equations (10) and (11) due to the excessive 
heating tends, initially, to override any effect of the upwelling. 
8y hour six there is only a .4 meter depth difference over that of S2, 
and essentially no temperature difference in the mixed layer. After 
that point, the upwelling begins to combine with the mixing to augment 
the cooling of the layer, and slowing its deepening. 
The net result after twenty-four hours is to produce 
a layer that 1s about as deep as that of S2, but that is as cool as 
the case of the evening upwelling. This represents a substantial 
energy difference in the runs. Again, a slight surface temperature 
Signature could be seen with this run in comparison with the standard. 
3. Gasenl II 
The next test 1s that of the interaction of mixing with a 
single cycle internal wave. B8y stepping the phase of an input sinusoid 
through the twenty-four hour progression of the model, a quasi two- 
dimensional motion can be prescribed. 
For a wave of the form 
wave amplitude ZA sin (wt + kx) 
the vertical motion would be of the form 
w(x,Z,t) @ Awcos(wt + kx) (21) 
as seen in Figure 6. 8y increasing the value of kx, the wave 1s 


made to travel in the -x direction. 
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The initial conditions used were those of cases S$] and $2. The 
vertical motion was introduced at hour one, and the amplitude was varied 
as t varied in equation (21), to give a wave induced vertical motion. 
This was then run twenty-four times, each time starting at a different 
phase of the wave, to test the interrelation of the wave input with the 
cyclical heat flux. For simplicity, only four major phase relation- 
ships are discussed here, as they will demonstrate the primary effects. 
Also, here, the two-dimensional field will be discussed for the first 
time. 

The wave form of the desired vertical motion does appear in all 
four of the displayed runs. The difference in each run is the phase of 
the vertical motion in relation to the heat flux (daily solar heating 
cycle). 

a. Shallow Initial Layer (Case IIIa) (Figures 27-33) 

In all four phase relationships presented, maximum downward 
motion at hours 1, 7, 13, and 19, the initial conditions are those of 
Sl. The magnitude of the peak upward and downward motions is 36 cm/hr 
at ten feet, as before. The wave is introduced at hour one, and pro- 
gresses through the system through hour twenty-four. 

(1) Peak Upward at Hour 7 (Figure 27). In this case the 
initial upward motion, hours one through twelve, the heating and the 
surface wind stress combine to give a slowly deepening layer, which, 
due to the upward advection of the layer and of cooler, below layer 
water, is almost .1°C cooler at its maximum temperature than SI. At 
the end of the unward cycle the layer nas deenened about two meters 


less than S1 at the same time, and has set up a slightly larger AT at 
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the entrainment zone. As a result, the following downward cycle pro- 
duces a smaller effect in combination with the wind stress in deepen- 
ing than that seen previously. Because of the initial upward motion 
effects on the temperature profile, the combined effect of wind stress 
mixing, surface cooling, and the downward motion produce a layer that, 
after twenty-four hours, is only two meters deeper than S1. Addition- 
ally, the downward effect of the motion on the temperature gradient 

is sufficient only to make up about 60 percent of the deviation in 
temperatures caused by the upwelling motion, as to that of Sl. 

The maximum 8T for this case is .245°C. It, also, 
produced the deepest layer of the four runs discussed here, over the 
twenty-four hour period. 

(2) Peak Upward at Hour 1 (Figure 28). In this case the 
upward motion was present during the first six hours of the heating 
cycle only. The downward motion cycle tgen progressed for twelve hours, 
followed by six hours of upward motion to complete the twenty-four 
hour sequence, 

During the first six hours the layer deepened to within 
less than one meter of the depth reached by hour six in SI, but the 
layer temperature is seen to be about .01°C cooler. The reversal of 
the motion was too late in the heating cycle to overcome that trend, 
and by hour nine, the time of the maximum surface temperature for Sl, 
this experiment produced a surface temperature that is .013°C cooler. 
The layer depth reached here is within half a meter of the S] depth 
for the same time. The difference is attributed to the initial lift- 


ing of the entrainment zone. 
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As the downward motion continues, along with the sur- 
face cooling, the layer deepens more quickly while the cooling of the 
layer slows as the effect of mixing is reduced. A maximum depth is 
reacned two hours after the downward motion has stopped, when the 
upward motion is of sufficient magnitude to overcome the wind stress 
mixing. The layer becomes shallower from its extreme depth of 15.56 
meters to 14,94 meters by hour twenty-four, with associated layer 
cooling due to surface heat loss, and increased entrainment due to 
the increased effect of the wind stress mixing. 

The maximum 8T for this run is .189°C. The final 
layer depth is within .3 meters of Sl for hour twenty-four, The 
variations of depth and temperature during the run are obviously 
much greater than for Sl. 

Peak Upward Motion at Hour 19 (Figure 29), The fact 
that the downwelling and the heating cycle are in phase, as in Case 
Ila, produces, by hour nine, a layer that is two nears deeper and 
.006°C warmer than S1 for this same time. The following upward motion 
is sufficient to halt layer deepening at hour sixteen in opposition 
to the mixing, but is not sufficient to cool the layer as much as for 
Sl. The final layer temperature is .035°C warmer than that in Sl, 
while the final layer depth is about one meter less, thougn at hour 
sixteen a depth of 14.2 meters was reached. The maximum 87 for tne 
mixed layer in this case was .181°C. 

(4) Peak Upward Motion at Hour 13 (Figure 30). This 
case produced the warmest layer temperature at hour nine, though the 


downward motion was only present during the first six hours of the 
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surface heating cycle. The initial downward motion deepened the layer, 
and thereby reduced mixing in this period, over the S$] case, allowing 
the surface layer to be heated to 14.988°C, an increase of .01°C over 
that in SI. The following upward motion was sufficient, hours eleven 
through fourteen, to hold the mixed layer at a relatively constant 
depth. This allowed a relatively constant, and strong wind stress 
entrainment to combine with Surface cooling to rapidly cool the layer. 
The result after twenty-four hours is to produce a layer of about the 
same depth as that of SI, but one that is slightly cooler. 

The 8T here was .237°C, a very noticeable difference 
from SI, though the final product did not vary much in layer temperature 
or layer depth, 

(5) Two Dimensional Data (Figure 31, 32 and 33). The two 
dimensional oceanic field was produced as previously discussed. The 
data covered here is that of the evolution of the ocean system at 
hours 6, 12, 18, and 24, 

By nour six the wave has progressed into the oceanic 
field, and its effect below the layer is dramatic, (Figure 31a). It 
has produced a mixed layer that has a mean depth of about 9 meters 
(initial h was 4.5 meters) witn a mean temperature of 14.957°C. A 
temperature structure, aS could be measured at the surface in x, 1s 
being formed. The present maximum difference in surface temperatures 
fer x) 1S" OTC , 

By nour twelve (Figure 31b) the effect of the wave is 
at its maximum, as seen in the amplitude of the contours below the 
layer. The layer is generally deeper than at hour six with a mean 
depth of about 11.5 meters. The deviation from the mean is much 
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greater for this time, however. The mean of the layer temperature 
has gone to 14,950°C, a slight decrease from hour six, but more 

Significant is the deviation of the surface temperatures in the x 
direction. There is now a very definite variation in x of .051°C. 

The layer has deepened to a mean depth of about 13 
meters by hour eighteen (Figure 3lc), with the effect of the wave 
beginning to be reduced in the system. The mean layer temperature 
nas dropped to 14,846°C, a cooling of almost .05°C in the six hours 
between illustrations. The surface temperature deviation has reached 
.091°C. Warm and cold temperature bands are becoming defined in x. 

The single wave has passed through the ocean by hour 
twenty-four, and only the remnant of its passage is seen (Figure 3la). 
In an ocean that is defined by SI conditions everywhere, the twenty- 
four hour product would be an ocean with a mixed layer depth constant 
at 14.65 meters, and mixed layer temperature of 14.767°C. With the 
introduction of a single diurnal wave cycle, a variable condition 
ocean has been produced. The mean depth of this ocean is 15.012 
meters, with local variations of up to +1.73 meters and -1.32 meters, 
with depths at small x being somewhat greater than depths at large x. 
The mean layer temperature is 14.764°C, but more importantly, the 
surface temperature variation is now ,095°C, in a clearly defined 
pattern. The pattern that exists is that of a narrow zone of rela- 
tively cool water centered near kx = 5, with a wider zone of relatively 
warmer water centered near kx = 18. Below the layer the isotherms are 
found to be parallel after the passage of this single wave cycle. 

b. Deep Initial Layer (Case IIIb) (Figures 34-40) 


The initial conditions for these runs were the same as 
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for S2. The wave cycle passed through the system as in the previous 
runs (IIIa). The wave form is harder to discern in these runs, but 
looking at the contours below the original layer (Figures 34-37) the 
same forms as the previous run are apparent. 

(1) Peak Upward Motion at Hour 7 (Figure 34). As in the 
previous deep initial layer cases, the shallowing due to the surface 
heating drives the model in the first six hours. The coincident up- 
ward motion during this time combines with the shallowing and mixing 
to produce a layer that is about 1.5 meters shallower, but at the 
Same temperature as in S2, by hour ten. Two noticeable effects during 
this period of upward motion are (i) a slight shallowing of the newly 
formed surface layer, and (ii) the upward transport of cooler below 
layer water, setting up a larger AT at the entrainment zone for this 
case over that of the standard. The downward motion then contributes 
to the deepening of the layer as does the mixing which must initially 
act over a relatively shallow depth. This entrains cool water, form- 
ing a quickly deepening and cooling surface mixed layer. As the down- 
ward motion decreases the cooling slows, and the increase in the depth 
of the layer is seen to be caused mostly by the vertical motion. Of 
all the phases examined this produced the deepest layer (30.34m) and 
has the second coolest layer at the end of the period (15.059°C) in 
relation to the other three cases discussed. 

The 8T for this case is .115°C. This is about nalf 
that for the shallow layer case discussed previously. The original 
deep interface can be followed on the contours (Figure 34b) over the 


twenty-four hour period. During the downward motion it, again, 1s 
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diffusing, while during the upward case there is a slight reinforce- 
ment of the interface. 

(2) Peak Upward Motion at Hour 1 (Figure 35). Again it 
1s seen that the shallowing mode of the mixed layer model, equations 
(10) and (11), controls the layer development to hour ten. The layer 
depth reached through the shallowing, upward, and downward motion 
combination by hour ten has reached a depth that is about .5 meters 
deeper than that in S2 for the same time. The temperature is the same 
as S2. Also, the downward motion makes entrainment of water into the 
mixed layer due to mixing more difficult, such that cooling after hour 
ten progresses much more slowly than in that of the previous case, and 
that of S2. The net result in the layer depth of this run is a layer 
that is about as deep as that of S2 after the twenty-four nour period. 
The model did not deepen more than S2 due to the dominant upward motion 
during the last six hours. The mixed layer temperature was sligntly 
warmer than that of S2 at the end of the run (.071°C). Also, it can 
been seen that the below layer water is somewhat warmer than in S2 for 
the end of the period, due to the combined effects of the heating and 
downward vertical motion timing. 

The 81 for this case is .092°C. This is less signi- 
ficant than in the previous run. 

(3) Peak Upward Motion at Hour 19 (Figure 36). Again, 
the shallowing mode of the model controls the initial ten hour layer 
temperature adjustment. The downward vertical motion, however, drives 
the layer two meters deeper during this time causing reduced mixing 


due to the surface wind stress. The upward motion keeps the layer 
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from deepening as quickly as before after hour twelve. Because of 
the below layer temperature adjustment by the downward motion and 
the initial deepening of the layer, the layer, also, does not cool 
as quickly. This case produced the shallowest final layer depth, 
and the second warmest layer at the end of the period, relative to 
the other three runs. 
The 8T for this case is .103°C, which is about the 

Same as the first case which was 180° out of phase with this run. 

(4) Peak Upward Motion at Hour 13 (Figure 37). The 
coincident shallowing, wind stress mixing, heat flux, and cycled 
vertical motion over the first ten hours produce a surface mixed 
layer that has the same temperature as S2, and almost the same layer 
depth (to within less than .5 meters). Below the layer, however, 
the upward motion has made relatively cooler water available by hour 
ten (in comparison with S2), and the layer, though not rapidly deep- 
ening due to the continued upward motion, cools quickly as large 
quantities of water are entrained over a relatively shallow depth, 
a trend which continues to the end of the period. The net result 1s 
a mixed layer that is 1.5 meters deeper than that of S2, but is also 
the coolest layer of the four runs. Because of the below layer temp- 
erature profiles which are set up by the upward motion for the mid- 
twelve hours, the continued run of this model over a second day gives 
continued rapid cooling of the layer due to increased entrainment of 
cooler, upward advected water, and to increased effects of mixing. 
A noticeable deepening does take place at the end of the period due 
to the downward motion, which accounts for the fact that this run 
ended deeper than S2. 
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The 8T for this run is .125°C. This is the largest 
of the four runs, but it is still only half of the maximum 8T seen for 
the shallow initial layer case. 

(5) Two-Dimensional Data (Figure 38, 39 and 40). The two- 
dimensional oceanic field is produced, here, as previously discussed, 
for the single wave cycle. The data covered is for hours 6, 12, 18, 
and 24, A common feature is the varying below layer temperature struc- 
ture set up because of the introduction of the initial deep layer. 
Because Of this, the changes in the thermocline below the initial deep 
interface best show the effect of the wave progression through the 
system. 

As in the shallow case, by hour six, the wave has 
progressed into the oceanic system producing some horizontal sea sur- 
face temperature structure. However, the effect on the surface 
temperature is almost negligible, in that the variability of the 
Surface temperature in x is only .001°C. This 1s due to the shallow- 
ing mode response of the model over the first six hours. Even with 
this, however, the depth structure has varied (Figure 38a) so that 
there has been a change from the constant structure that would be 
seen in S2 of 11.767 meters for the layer depth, to a layer of mean 
depth of 11.769 meters J .35 meters, indicating that an adjustment 
due to the wave is taking place. 

By hour twelve much more temperature and layer (h(x)) 
variation is beginning to take shape. At hour twelve the case without 
the vertical motion has a layer depth of 15.037 meters, and a layer 
temperature of 15.159°C everywhere. In this case, where the maximum 
effect of the wave has been reached, the layer depth has a mean of 
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15.231 meters with a deviation of t 2.5 meters. The variability in 
the surface temperature in x has reached .01°C. Though the variabil- 
ity in the surface temperature is not very large, a pattern is again 
beginning to emerge in the x direction. At small x cooler water is 
forming, while at large x a band of warmer water is emerging. 

The layer has continued to deepen to hour eighteen. 
The layer depth has increased to a mean of 19.46 m, * 3.4 m, and the 
variation in the surface temperature structure in x has reached .03°C. 
The warm-cold pattern seen at hour six is repeated and strengthened. 
Though the temperature difference between warm and cool water is not 
as great as in the shallow case, the banding of the patches is much 
more distinct here (Figure 38, Figure 31). 

At the end of the period the wave has, again, been 
removed from the system, and only the remnant of its passage is seen. 
The final layer depth is deeper than in the previous hours, and much 
deeper than in the shallow case. The variability of h in x is also 
much greater here than in the shallow case. The surface temperature 
patchiness is well-defined, but the variability of the temperature in 
x is omly .035°C. 

Comparing the two cases, shallow and deep initial 
layers, it is seen that for the same amplitude wave the model produces 
two distinct ocean types. For the shallow case an ocean is formed 
that has a relatively shallow layer with a large temperature varia- 
tion in the x direction, but with that temperature variation primarily 
in wide bands. For the deep case a deep surface layer is formed with 
smaller variability in the surface temperature. However, the temp-. 
erature structure in x iS much better defined. Also in the deep case, 
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the patches of warm and cold water are much closer to each other (in x). 
This could suggest a mechanism for generating or sharpening SST fronts 
through surface heating, and internal solitary wave passage. 

4. Case IV 

It was noted in Case III that after twenty-four hours the wave 
form had passed through the ocean system, and had left behind only the 
residue of its effect on the system, over and above the standard effects. 
This can be compared to a solitary wave passing through an ocean system. 
While this has theoretical advantages, it is not a very realistic prob- 
lem, in that wave forms seldom take on this sort of motion. One pos- 
sible exception may be solitons generated at extreme tidal stages 
(Apel, 1982). 

The next progression, therefore is to look at a continuous wave 
system that passes through the ocean with constant form, until the run 
is terminated. This is done by simply modifying the mathematics of 
the previous case to include initialization of the temperature and 
salinity profiles with the wave form according to the following. 


We know that 


w(x ,Z,t) 


AAzcos(wt + kx) 


and 


vy = -@ 3h 
dt 


or 
Sh = - f war 
therefore 


Sh —=- -AAzsin(wt + kx) 
wD 
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For initial displacement 


8h(x,z,0) = Sh 


therefore 


Sho = —AAzsin(kx) 
wD (22) 


Applying equation (22) to form the initialized layer depth, and com- 
bining this with equation (16) to initialize the temperature profile, 
a continuous wave field can now be passed through the model ocean 
(Figure 7). 

From this point on the discussion will be confined to the deep 
initial layer case, S2 standard, and variations thereof. This is done 
because all the primary effects discussed to this point appear in both 
the shallow and deep cases, as do the effects to be discussed in vary- 
ing degrees, and it is not necessary to continuously reiterate this 
discussion, 

a. Diurnal Wave (Case IVa) (Figures 41, 42 and 43) 

Because of the overriding effect of the shallowing mode of 
the model in the first six hours of the runs, the results produced 
here were exactly the same as those discussed in Case IIIb, for the 
phase relationships. This was of course expected from the theory. 

(1) Two-Dimensional Data (Figure 44, 45 and 46). Again, 
because of the dominance of the shallowing mode, the hour 6, 12, 18, 
and 24 two-dimensional fields produced here are the same as those of 
Case IIIb for the area of the entrainment zone, and the mixed layer. 
Below the layer, however, the wave form 1S now seen to progress through 


the system in a continuous mode. This is best seen in the contours 
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in which the amplitude of the system remains a constant, but the wave 
form is shifted with time. 

It can be seen ete the patch of warmer waters forms 
(Figure 44) and takes a phase that barely leads the peak of the wave 
system. The entrainment zone ( h(x) ) now takes the form of the 
passing wave system, as would be expected, and yet was not obvious in 
the previous case. The divergence and convergence of below layer 
fields 1S also obvious here, where previously it was not particularly 
obvious. 

b. Semi-diurnal Wave (Case IVb) (Figures 47-50) 

A dominant internal wave system of the temperate oceans 
1s the semi-diurnal tidal cycle. South of 30° N there is also a 
strong diurnal tidal cycle. A semi-diurnal tidal field can be pro- 
duced for testing in this model by simply adjusting the frequency of 
the prescribed w for equation (21), and setting the initial profiles 
with the appropriate frequency of equation (22) and (17). 

The amplitude of the motion remains the same as the 
diurnal cycle case. Because the cycle of wave mixing interactions 
repeats itself after a 2m phase shift, which occurs after 12 hours 
in this case, only the cases of six and twelve hour phase shifts 
will be discussed. The 18 and 24 hour phase shifts are merely repeats. 

(1) Six Hour Shift (Figure 47a and b). Figure 8 shows 
the relationship of the heating cycle to the wave motion i this case, 
The vertical motion is seen to immediately combine with the shal lowing 
mode of the model through the initialization to give a somewhat 
shallower initial layer. The shallowing mode of the model holds the 
layer temperature the same as in S2 to hour six, while the vertical 
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motion varies the layer depth, and alters the below layer temperature 
Structure. B8y hour six a half-cycle of the wave has passed through 
the model point, and a layer has been produced that is almost one 
meter deeper than S2 at this time, but which has the same layer temp- 
erature. Continued but decreasing downward motion combines with the 
downward heat flux until hour nine to slowly deepen the layer, but 

the layer continues to warm. This process continues through hour ten, 
though upward motion has begun because the heating and entrainment of 
relatively warm water are greater than the effect of the resistance to 
deepening due to the surface heat flux, and the small upward motion. 
By hour eleven the upward motion dominates, and the layer shallows 
while it cools slightly due to the increased entrainment. At hour 
twelve the layer is as shallow as it will ever get, and the first wave 
has passed. Now decreased upward motion followed by downward motion, 
combined with surface cooling and mixing through a progressively cooler 
surface mixed layer quickly increases the layer depth to a maximum of 
26.7 meters at hour twenty-one. This is deeper than S2 gets in the 
entire twenty-four hour period. Yet, the layer is not as cool as in 
the S2 case. 8y hour twenty-one the upward motion is again dominant 
and the layer shallows. 

(2) Twelve Hour Shift (Figure 47c and d). Figure 9 shows 
the relationship of the heating cycle with the wave in this case. The 
initialization combines with the vertical motion and shallowing mode 
to produce an initial layer that is deeper than that in S2, but soon 
becomes shallower, as the new surface layer is established. By hour 


six a layer of the same temperature as S2, but ./2 meters shallower 
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has been formed. The net result of the cyclical interactions in this 
case is a layer that is 1.5 meters deeper than in S2, but which is at 
about the same temperature. The difference in the layer temperature 
between this run, and the previous is .004°C. However, the layer 
depth is 2.2 m different. 

(3) Two-Dimensional Data (Figure 48, 49 and 50). At hour 
Six a very discernable depth structure has been set up at the bottom 
of the mixed layer which reflects the structure of the passing wave. 
Above the entrainment zone the temperature structure, as produced by 
the shallowing mode, is homogeneous. Below the layer the structure is 
highly varied in x reflecting the combined effects of the vertical 
motion to this time. The mean layer depth is 11.75 meters, while that 
of S2 is 11.77 m, showing the domination of the shallowing to this point. 
A variation of h in x of * .75 meters is already exhibited. 

By hour twelve the mixed layer depth structure in x, 
and the layer temperature variability in x have evolved considerably 
over that of hour six. At the entrainment zone the structure is approx- 
imately sinusoidal, as expected. Over the last six hours a large 
amount of water from below the layer has been entrained while surface 
heating has occurred, and the mixed layer has heated and cooled 
accordingly to begin to set up the pattern of warm and cold bands of 
temperature in x. The mean depth of the layer is now 15.2] meters 
* 2.50 meters, while that for S2 is 15.04 meters. The temperature 
variation in x as seen at the surface is now .012°C. 

The structure of temperature and layer depth profile 


has been further accented by hour eighteen. The layer has continued 


77 








to deepen, and the fluctuation in x increased. The banding of the 
in layer temperature pattern in x has become more defined. 

After twenty-four hours a layer has been produced 
that has a depth of 26.14 m, * 3.3 m. The surface temperature pattern 
in x 1S quite clearly defined, though the variation in temperature is 
only .015°C. 

5. Case V (Superimposed Waves) (Figure 51, 52 and 53) 

Since oceanic internal waves ae almost never simple sinusoids, 
one more progression is considered. Spectral analysis of the MILE 
data (Figure 10) at station P (55° N) shows the dominant wave field to 
be the semi-diurnal tide superimposed by, at least, the 6-7 hour har- 
monic of the tidal wave. Internal wave theory indicates that the 
period of the internal waves must be less than or equal to the inertial 
period. North of approximately 30° N, no diurnal internal waves may 
propagate, and therefore, they are not seen in the MILE data. South 
of this, however, diurnal and semi-diurnal tidal period waves may 
Propagate. Therefore, since previous discussions have been concerned 
with diurnal and semi-diurnal waves, and there is justification for 
their presence, these waves will be superimposed to form a slightly 
more complicated wave system. 

The combined wave is directed through the two-dimensional 
ocean which has the initial conditions of S2, from right to left 
(large kx to small kx). The amplitude of the motion is such that 
both waves peak at the same amplitude at their peak motions, though 
one has half the period of the other. Both waves start with an 


amplitude of vertical motion in phase at time one. 
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The superimposed waves reflect all of the interactions discussed 
in the previous sections. Further, this simple addition of two wave- 
forms has developed several new features. The new features can be 
traced in their development from hour one through hour twenty-four, 
where they become most evident. Since this is the case, the discussion 
will be limited to the case of the two-dimensional ocean at hour twenty- 
four. 

By this time a complex structure has emerged for the two-dimen- 
sional ocean, both in temperature and in profile of the layer depth. 

An ocean has been formed that has a mean layer depth of 26.78 m, but 
this depth has a variation that is not symmetrical from small to large x. 
At small x the variation from the mean depth is seen to be a maximum of 
-5.45 m, while at large x the layer depth varies to a depth deviation 
from the mean of +7.83 m. An ocean has been formed that has a deep 
layer at one end, and a shallow layer at the other. Further, a devia- 
tion of the surface temperature in the x direction of .042°C is seen. 
The pattern of the warm and cold bands again is seen, but this also is 
no longer symmetrical. The temperature pattern, as seen at the surface, 
is a cold band beginning at small kx, progressing through a weak gradi- 
ent to a warm band centered near kx = 15. A strong gradient is then 
encountered leading to a cold band centered about kx = 21. Through the 
introduction of a simple two wave system, an extremely diverse ocean 

has been produced. 

Shifting the phase relationships of the two internal waves, and 
running the model again, resulted in similar patterns. Positions of 
the temperature maxima and minima in the layer, and the position of 
the depth variations shift, but their relative positions are maintained. 
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The temperature variation in the x direction as seen at the surface 


does not change with the phase changes. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


The analysis of the mixing response to the prescribed vertical 
motion shows that vertical advection can be successfully introduced 
into this one-dimensional mixed layer model without numerical diffi- 
culties. This is important in that it indicates that a single station 
forecasting model can be altered to include the variability of long 
internal wave motion. This step is the first of a progressive devel- 
opment of a single station forecasting model that may in the future 
include other local advective effects. 

Using an ensemble of one-dimensional stations, a two-dimensional 
mixing response to prescribed internal wave motion was simulated. 
Although the initial temperature field (without wave) and the surface 
boundary conditions were both horizontally homogeneous, a significant 
wave-induced horizontal structure was created in the T(x,z) field. Both 
linear and non-linear effects have been demonstrated in thé simple 
single wave and superimposed wave cases. Horizontal variability has 
been shown, both in mixed layer depth and in mixed layer temperature 
structure. The mixed layer temperature develops a corrollary structure 
of alternating warm and cold patches. In the case of the superimposed 
waves of different frequencies, this patchiness pattern can be complex, 
and can qaive rise to horizontal gradients which could be frontogenetic. 

There is an important relationship between the phase of internal 
motion and the cyclical diurnal heat flux boundary condition. This is 
clearly demonstrated in the numerical experiments in which the phase 
of the input pulse or wave system has been varied. 
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In the case of superimposed diurnal and semi-diurnal waves, the 
horizontal variation in temperature in the mixed layer was about .04°C 
after only a single diurnal cycle. In the same period, the mixed layer 
depth varied by about 13 m. The wave dispersion relation for a two layer 
ocean at 25° N gives a horizontal wave length scale of from 10 to 100 Km, 
depending upon stratification. This wave-mixing interaction may, there- 
fore, be an important mechanism for the generation of sea surface 
patchiness as found by Butler (1981). 

Further study must be done in this area. Clearly, more complicated 
waveforms will produce more complex responses. The model should be 
tested with salinity changes, including evaporation and precipitation. 
Tests should be conducted on data sets such as those of MILE. Compar- 
ison of the mixed layer model results using HILE boundary conditions, 
with and without vertical advection, as observed in spectrum analysis, 
should further demonstrate the importance of addition of vertical 
motion to mixed layer models. 

Finally, the wave systems should be tested in an interactive mode 
with the mixing of the model. Energy exchanges, as prescribed by wave 
theory and the studies of Linden (1975), and Kantha et al. (1977) 
should be allowed, and the relative significance of the different 
mechanism should be evaluated. Internal vertical motion is an impor- 
tant dynamic variable, and should be treated as such in the future of 


mixed layer modelling. 
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1 i? 09s2 15.112 34.000 0.50 100.8 0.0 0.0 
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Same as figure (1), except 
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DAY HR 4 t s TAU 2NET WDD OH 
lol $1.3 15.001 34.000 0.50 = 55.7 0.03600 191.500 
io 16.6 15.009 34.000 0.50 -201.5 0.03600 65.705 i 
Oe a ee 15.024 34.000 0.50 -326.8 0.03600 38.9846 | 
1 $4 13.84 18.028 34.000 0.50 -423.0 0.03600 47.351 4 
2 5-12.48 15.076 34.000 0.50 -483.5 0.03600 43.922 
a ane ee 15.105 34.000 0.50 +504.2 0.03600 42.862 Te) n = 
9) Se 15.133 34.000 0.50 -483.7 0.03600 44.985 10 fh aS 3 
bo) 8 1355 15.155 34.000 0.50 -423.5 0.03600 47.730 ‘A x 
173 14.4 15.170 34.000 0.50 -327.5 0.03600 §1.043 a ~ 
lL 10 15.6 ~- 15.176 34.000 0.50 -202.4 0.03600 55.035 | 
eee Ge 15.175 34.000 0.50 = $6.6 0.03600 $9.415 
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1 ars 6 15.157 34.900 0.50 100.8 0.03600 69.859 | 
bee 2be2 15.149 34.000 0.50 100.8 0.03600 75.211 
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Lie 26.5 15.130 34.000 0.50 100.8 0.03600 91.896 | 
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= 3 
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OaY HR 4A T 3 TAU QNET wDD DE 
ed 6.0 14.956 34.000 0.50 - 3.5 -.03600 - 22.140 
1 <2 6.9 14.937 34.000 0.50 -100.8 -.03600 - 25.258 
1 3 7.3 14.931 34.000 0.50 -184.3 -.03600 - 26.750 
1 4 7.6 14.934 34.000 0.50 -248.4 -.03600 - 27.962 
1 5 7.8 14.941 34.000 0.50 -288.7 -.03600 - 28,725 
l 6 7.9 14.949 34.000 0.50 =302.5 -.03600 - 28.917 
eS? 7.9 14.957 34.000 0.50 -288.9 -.03600 - 29.021 
1 & 7.9 14.961 34.000 0.50 -248.7 =.03600 - 29.089 
1 & 8.0 14.959 34.000 0.50 -laa.7 -.03600 = 29.169 
1Lo 8.0 14.949 34.000 9.50 -101.3 -.03600 - 29.288 
fenti 8.0 14.930 34.000 0.50 = San -.03600 = 29.467 
i-vic a.1 14.898 34.000 0.50 100.2 -.03600 - 29.9817 
1°23 3.3 14.869 34.000 0.50 100.8 -.03600 - 30.262 
1 14 8.4 14.342 34.000 3.50 100.3 -.03600 - 30.793 
1 15S 3.6 14.915 34.000 0.50 100.3 -.03600 - - 31.417 
1 16 8.7? 14.789 34.000 9.50 100.38 -.03600 - 31.994 
ee a7 14.763 34.000 9.50 100.8 -.03600 - 31.992 
1 18 a. 14.735 34.000 0.50 100.8 -.03600 = 31.779 
Pais 9.7 14.709 34.000 0.50 100.8 =.03600 = 31-733 ~ 
Los $.? 14.683 34.000 0.50 100.8 -.03600 - 31.750 > 
lo2 $52 14.657 34.000 9.50 100.8 -.03600 = 31.832 — 
lore2 8.6 14.631 34.000 0.50 100.3 -.03600 = 31.625 
1 23 3.6 14.605 34.000 9.50 100.8 -.03600 - 31.475 x= 
I 24 3.6 14.573 34.000 0.50 100.8 -.03600 - 31.394 ms 
- 1 8.5 14.556 34.900 9.50 - 3.5 -.03600 - 31.297 LJ 
2 2 8.5 14.547 34.000 0.50 -100.8 -.03600 - 31.155 =) 
ee a5 14.546 34.900 0.50 -194.3 -.03600 - 30.956 
eer | 9.3 14.551 34.900 9.50 -248.4 -.03600 - 30.709 
aces 3.3 14.559 34.000 0.50 -288.7 -.93600 - 30.435 
2 6 9.2 14.554 34.000 0.50 -302.5 -.03600 - 30.162 
22 3.2 14.562 34.000 9.50 -288.9 -.03600 - 29.951 
293 8.1 14.567 34.000 0.50 -248.7 -.03600 - 29.802 
2 9 9.1 14.565 34.000 0.50 -184,7 -.03600 - 29.745 
2. 10 aol 14.555 34.000 0.50 -101.3 -.03600 - 29.308 
Pa 3.2 14.535 34.000 0.50 = 4.1 -.03600 - 30.018 
2 12 8.3 14.477 34.000 0.50 100.2 =.03600 - 30.407 
Pe | 3.5 14.449 34.000 0.50 100.8 -.03600 - 31.019 
2 14 9.7 14.422 34.000 0.50 100.8 =.03600 - 31.729 
2 #15 8.7 14.396 34.000 0.50 106.3 -.03600 - 31.972 
2 16 a.5 14.369 34.000 0.50 100.8 =.03600 = 31.256 
va Ne 8.3 14.341 34.000 9.50 100.8 -.03600 - 30.583 
2 #18 8.2 14.304 34.000 0.50 100.8 -.93600 - 29.953 
2 #19 8.0 13.274 34.000 0.50 100.8 =.03600 - 29.370 
2 20 7.9 14.243 34.000 2.50 100.8  =.03600 - 28.928 
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2 23 4 14.146 34.000 9.50 100.8 -.03600 - 27.302 
2 24 73 14.095 34.000 0.50 100.8 -.03600 = 26.852 
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wie 27.6 15.166 34.000 0.50 100.8 0.02546 69.477 
22 28.9 15.162 34.000 0.50 100.8 3.91800 £1.644 
23 30.0 15.157 34.600 0.50 106.8 0.00932 27.833 


nal i 

OAY HR 8 t $s TAU QNET woo OK i i 
ue get 49.5 15.001 34.000 0.50 - 55.7 0.0 0.0 - : 
1 2 To. 2 15.009 34.000 0.50 -201.5 0.0 0.0 5 
1 3 ic 2 15.024 34.000 0.50 =326.8 0.0 0.0 i 8 
Ll 4 13.0 15.04? 34.000 0.50 +-423.0 0.0 0.0 5 rn 
1 5 a. 1 15.075 34.000 0.50 -483.5 0.0 0.0 "| 1 
Ll 6 11.8 15.105 34.900 0.50 ~-504.2 0.0 0.0 10 2 
aa Tm. 1 15.132 34.000 0.50 -483.? 0.0 0.0 fi — 
1a 12.5 15.154 34.006 0.50 423.5 0.0 0.0 t wn 

1 93 13.0 15.168 34.000 0.50 -327.5 0.0 0.0 | a 

i toe 3.5 15.175 34.000 0.50 -202.4 0.0 0.0 § 

1 wire 4.2 15.172 34.000 0.50 - 56.6 G.0 0.0 i 

1522 15.6 15.160 34.000 0.50 99.9 0.0 0.0 

1 13 16.0 15.149 34.000 0.50 100.8 0.00932 14.869 

t °2 722 15.139 34.000 0.50 100.8 0.01800 30.739 

1 &S 18.6 15.130 34.c00 9.50 100.8 0.02546 46.756 

ieee. 0 15.122 34.000 0.50 100.8 0.03118 61.334 20 

Ne) se a 15.114 34.000 0.50 100.8 0.03477 73.399 : 
Tete 2 2, 2 15.107 342000 e.50 100.a 0.03600 91.528 

1° oie. 24.7 15.101 34.900 0.50 100.8 0.03477 84.247 

I zo" 26.2 15.096 34.000 0.50 100.8 0.03118 80.442 

Leet 2? 15.090 34.500 0.50 160.8 0.02546 69.511 = 

1 22 28.9 15.986 34.000 ¢.50 100.8 0.01800 51.613 =< 

1 =: 30.0 15.081 34.00¢ 9.50 100.8 0.06932 27.799 

1 24 30.8 15.077 34.900 0.50 100.8 0.0 0.0 B 

fe 

ees el els oo 734.000 0.50 = 55.7 0.0 0.0 af q 
2a 21.0 15.085 34.000 0.50 -201.5 0.0 0.0 lJ 

23 15.4 15.101 34.000 0.50 -326.8 0.0 0.0 Qi 

eee Tesh 2 1S.3Zz7> 34.000 G.50 ~-423.0 0.0 9.0 cz 

2 5 12.0 15.152 34.600 0.50 -483.5 6.6 0.0 

2 6 Co. 7 15.181 34.000 0.50 -+504.2 0.0 0.0 

2 12.0 15.208 34.000 0.50 -483.? 0.0 0.0 

2 9 12.4 15.231 34.0900 C.50 -423.5 o.¢ 0.0 

2 9 12.9 15.245 34.000 0.50 -327.5 0.0 0.0 t= 

ee te ee Oe | 15-251 382c6C 0.50 +-202.4 0.0 0.0 40 

ye eg a 15.248 34.000 9.50 - 56.6 9.0 9.0 

2s 14.9 15.236 34.000 0.50 39.9 0.0 0.0 E 

Zoids, 15:9 15.225 34.000 0.50 100.8 0.00932 14.752 

a iae 17.0 15.215 34.000 ¢.50 106.8 0.01800 30.358 

Zurtse 19.3 15.206 34.c00 °.50 100.8 0.02546 46.032 

2 #16 19.8 15.198 34.000 4.50 100.6 0.03118 60.690 — 

See 23,3 15.190 34.000 0.50 100.8 G20sa 07 Te. 7194 o 
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sz, AYER Teer. AYER JEPTH DE: 
5.3009 49.5152 
.5.1690 19.1073 . | 
55.0280 16.1826 
or 15.0470 13.0180 ( Cc ) | 
| 1$.0750 12.0605 | 
ass | 15.1040 11.7679 | 
-. »§.1310 12.1436 
tei 5.1530 .2.4a81 | 
15.1690 32.9925 
5.1740 13.5268 
o ‘$4719 14,2001 
S | "5.1590 15.0372 
. 15.1480 26.0323 } 
2! ‘S280 17.2318 | | 
—4 15.1290 18.5022 a 
25.1219 19.9813 | 
5.1140 21.4770 | 
15.1370 23.957) 
=) 5.1039 24.5510 
=) 5.0950 26.2263 | 
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= .5.5350 28.9314 | 
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, Figure 19. Same as figure (1), except 
for Case IIb wath peak at hour 18. | | 
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0.50 
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9.50 
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0.50 
0.50 
G.50 
0.50 
0.50 
0.50 
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- $5.7 
-201.5 
326.8 
~423.0 
“483.5 
-504.2 
-483.7 
-423.5 
327.5 
7202.4 
- $6.6 
99.9 
100.28 
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100.8 
100.8 
100.8 
100.8 
106.8 
100.8 
10c.8 
100.8 
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1¢c0.38 


LG & 
-201.5 
-326.8 
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Same as figure (1), except 


OAY HR tT $ 
2 S$i.3 15.001 34.000 
1 2 18.5 15.009 34.000 
1 3 16.5 15.024 34.000 
1 o4 13.2 25.048 34.000 
L & Yami 19.076 14.500 
1 6 11.8 °%5.105 34.000 
, 2 iat 18-132 442000 
* @ WES 15.1546 34.000 
1 9 13.0 15.169 34.000 
1 20 13.6 15.175 34.000 
L GY eee «86 1S sa72 «6034000 
1 12 15.1 15.160 34.000 
1 13 15.3 15.149 34.000 
1 14 16.8 15.139 34.000 
L 1S 17.6 15.130 34.000 
1 16 18.5 15.121 34.000 
Tota ses )=— 15-213 «34.060 
176 92622 18.105 34.000 
te F203 086 1S2098 34.500 
1 20 22.5 15.091 34.000 
L 22 24.0 15.985 34.000 
S220 25-7 i55079 346.000 
1 23 27.6 38.074 34.000 
1 28 29.5 15.069 34.000 
2 Lb 32.0 15.07¢ 34.000 
z 2 21.7 15.077 34.600 
2 3 15.€ 5.093 34.200 
2 § 413.3 15.116 34.000 
2 S$ 12.1 15.144 34.000 
26 22 15.173 34.000 
2 7 12.0 25.201 34.006 
Z 8 dz.4 82223 34.000 
2 3 12.9 18.236 34.000 
2.19 «#413.5 15.244 34.000 
2 12 14.2 15.241 34.000 
Sie 5.0 155029 34.¢00 
2 13 18.8 15.216 34.000 
2. 14 16.7 15.208 -34.000 
2 15 17.5 15.19% 34.660 
2. 16 18.3 15.189 34.000 
2 7 e238 152001 34.000 
2 WH 2651 15.173 14.000 
2 le 21-2 15.166 34.690 
3 3G 22.5 15.189 34.900 
2 @ 24.1 18.153 34.000 
+ 22 25.3 15.147 34.000 
252392727 15.182 34.600 
fa2@ 2@se i$223° 34/000 
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S 
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S 
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1§. 9600 $1.2673 
15.9096 19.%99i1 
15.0240 16.4907 
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: 1$.0750 12.1330 
- 1§. 1940 11.7510 
aa 15.1330 12.1393 
15.1589 12.4952 
15.1630 13.0005 
25.1740 13.5514 
S 15.1720 14.2418 
=) 1$. 1660 15.0aN4 
° 75.1490 15.3397 
2 15.3399 16.7917 
>4 25.1290 17.4410 
15.1210 16.4985 
| 1§.3120 {9.3379 
15.1350 20.1931 
2 15.0970 21.2566 
3 15.0910 22.5307 
F 15.9840 24.0162 
=} 1$.9799 25.7130 
=! 15.2730 27.5654 
= 15.0680 29.5215 
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Figure 20. 
= 
3 | 
*2 
an” 
Du 


bu. UU 


12.90 


12.50 


13.00 


-0 


i*.00 
(RELATIVE TO 


Ou 


175.207 


56.764 
41.448 
23.522 
11.257 


veooooeoeocoeooncdces 


eococnaoocaeoo0oeeco 


0.0 
19.714 
40.192 
60.364 
78.928 
94.207 


104.386 


106.008 


66.710 
39.216 
23.782 
11.182 
a.¢ 


eoocooensoocen 
oooooqnooecoecono 


0.0 
19.677 
410.216 
60.511 
79.246 
94.597 


104.837 


14.50 








a ce 


Z 


15.00 


FURST) 


Ss 


| 
| 


DEPTH (M) 





i$.50 





16.00 





i 
| 


67, | 





TIME CHR) 


(b) 


16.50 





DAY 4YR a3 $ TAU NET @00 Ou 
lot $0.4 18.0021 34.000 0.50 ~- $5.7 0.01800 39.930 : 
1 2 18.5 18.009 34.000 0.50 -202.5 0.02546 46.948 
1 3 16.6 1$.024 34.000 0.50 +326.8 0.03118 5$0.926 tt 
Lb 4 1354 19.048 34.000 0.5C +-+423.0 0.03477 45.750 
1s 12.4 18.076 34.000 0.50 -483.5 0.03600 43.926 
i 6 ze4 1S.195 34.000 0.50 -S04.2 0.03477 41.376 10 
9 12.6 1S8.133 34.000 0.50 -483,7 0.03118 38.787 
1 3 i | 18.155 34.000 5.50 -423.5 0.02546 33.483 
1 9 baat 18.169 34.000 0.5C -327.5 0.01800 25.150 
1 10 14.9 15.176 34.000 0.50 -202.4 0.00932 13.952 
1 ll 15.6 15.174 34.000 0.50 - $6.6 0.0 0.0 
Lb 12 16.5 18.163 34.000 0.50 99.9 o.6 0.0 
1 123 17.3 15.184 34.000 9.50 100.8 0.0 0.0 
1 la 38.1 18.145 34.000 0.50 106-487) 0.0 0.90 
1 1S 18.9 15.136 34.000 0.50 100.0 90.0 2.0 
1 16 19.7 18.128 34.000 ¢.S$0 100.8 0.0 0.0 20 
lo o17 20.5 15.121 34.000 0.50 100.8 9.0 0.0 
1 18 21.4 15.114 34.000 0.5 100.8 0.0 0.0 
1 le #222 18.107 34.000 0.50 100.8 0.0 0.0 
lL 20 23.0 18.100 34.000 0.50 100.8 92.0 9.0 
ae 23.9 18.094 34.000 0.50 100.8 0.0 0.0 - 
1 22 24.7 15.089 34.000 0.50 100.3 9.0 0.0 = | 
12d 2s.6 18.08} 34.000 0.50 100.8 0.6 2.0 ~, Z 
1 28 26.7 25.078 34..6¢€ 0.50 100.3 0.00932 24.743 Eff 
p= 
2 2 27.6 15.078 34.000 0.50 - $5.7 0.01800 49.312 a? 
22 2.6 1$.985 34.900 3.56 -201.5 0.02546 $6.700 uJ 
2 3 18.5 1$.102 34.900 0.50 -326.8 0.03118 47.578 a 
oe 13.6 18.124 34.000 0.50 -423.9 0.03477 46.411 
cs t2.3 1$.152 34.000 0.50 -483.5 0.03600 43.595 
2 6 20k 15.182 34.000 9.50 +504.2 0.03477 41.266 
23 228 15.210 34.000 9.50 -483.7 9.03118 38.240 L. 
2 3 Mas 2 18.232 34.000 0.50 +423.5 0.02546 332220 
225 14.0 18.247 34.900 0.5¢€ -327.5 0.01806 24.907 
2 #10 14.7 18.253 34.000 9.50 +202.4 0.00932 13.616 40 
2. Ui1) Se 1$.251 34.000 0.50 + $6.6 0.0 0.0 
2 12 “3622 18.241 34.006 2.50 99.3 0.0 0.0 
2003: el 15.231 34.000 0.50 Hod.8 0.0 0.0 
2) 214 (9 15.222 34.000 0.50 100.8 0.0 0.0 
2 1S 168.7? 15.213 34.000 0.50 100.3 90.0 0.0 
2 hfe 19.5 15.205 34.000 0.59 100.8 0.0 9.0 
2 by 20.3 415.197 34,000 0.50 100.8 0.0 0.0 
2 te 21.2 18.190 38.000 0.50 100.8 90.0 9.0 
2 Yee 22.0 15.283 34.000 9.50 100.8 0.0 2.0 Si 
2 w@@ 22,8 15.177 34.000 0.50 160. p 0.9 9.0 
2. 2) 25,7 15.171 34.000 2.50 100.8 90.0 0.0 
Pi 24.6 18.165 34.000 0.50 100.8 90.0 0.0 
2° 23°" 25.5 18.159 34.000 0.50 100.8 0.9 0.0 
a 2 26.7 18.154 34.000 0.50 100.8 0.00932 24.7 
i (a) 
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*$.347 13. 1868 
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~~ 5.:320 ‘2.6361 
5.1540 43.3217 
§. 830 4.0968 
{Soa30 4.938? 
a 5.1730 5.6259 
5S 1§.1530 26.4616 
. *5..533 7.2791 
bd S. 1440 8.0945 
7% §.1360 . 9. 3090 
.5..289 29.7240 
.$. 1206 23.5437 
:$.9 730 21.3563 
4 5.1360 22.1925 
S 5.1500 23.9234 
a $.2900 23.6532 
ie $.0990 24.6984 
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=o 
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= Figure 21. Same as figure (1), except ! 
) for Case [Ib with peak at hour 35. 
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OAY HR 4 T 3 TAU QNuET woo OH 
! 
ff 6.0 14.956 34.000 0.50 - 3.5 -.03477 -21.398 
1 2 6.9 14.937 34.000 0.50 -100.8 -.03118 -21.939 
1° °3 7.4 14.932 34.000 0.50 -184.3 -.02546 “19.111 
bo <4 7.9 14.934 34.000 0.50 -248.4 -.01800 -14.295 
1 5 8.2 14.943 34.000 0.50 -288,7 -.00932 - 7.661 
l 6 8.5 14.3950 34.000 0.50 -302.5 0.0 0.0 
YY 8.9 14.960 34.000 0.50 -288.9 0.0 0.0 
lea 9.2 14.966 34.000 0.50 -248.7 0.0 0.0 
1 9 9.5 14.967 34.000 0.50 -184.? 0.0 0.0 
1 10 9.8 14.962 34.000 0.50 -101.3 0.0 0.0 
Lee loon? 14.949 34.000 0.50 - 4.1 0.0 0.0 
l 12 10.4 14.927 34.000 0.50 100.2 0.0 0.0 
1 13 10.8 14.308 34.000 0.50 100.8 0.0 0.0 
el ey Staal 14.889 34.000 3.50 100.a 2.0 0.0 
ik YS llee 14.871 34.000 0.50 100.8 0.0 0.0 
Le Greek) 27 14.854 34.000 0.50 100.38 0.0 0.0 
lg 12.0 14.838 34.000 0.50 100.8 0.0 0.0 
Y ie) 12.3 14.821 34.000 0.50 100.8 0.0 0.0 
1 #19 22.5 14.806 34.000 0.50 100.8 -.00932 -11.687 
a2 Oe 2. S 14.791 34.000 0.50 100.8 -.01800 -22.791 
sh 2h) 2.5 14.776 34.000 0.50 100.8 -.02546 -32.299 = 13.8491 
zie 2.4 14.761 34.000 0.50 100.8 -.03118 -39.426 = 
leezjemi2. 3 14.746 34.000 0.50 100.3 -.03477 -43.687 —— 
1 24 12.2 14.727 34.900 0.50 100.8 -.03600 -44.903 oa 
2 1 12.2 14.719 38.000 0.50 = 3.5 -.03477 -43.006 a 
22 2rerk: 14.717 34.000 0.50 -100.8 -.03118 -38.216 ul 13.4981 
za 12.0 14.721 34.900 0.50 -184.3 -.02546 -30.974 OQ 
216 12.0 14.730 34.000 0.50 -248.4 -.01800 ~21.820 
2 = § PEGA 14.742 34.000 0.50 -288.7 -.00932 Sele sli 
2 6 eae 14.749 34.000 0.50 -302.5 0.0 0.9 
SG 12.4 14.761 34.000 0.50 -288.9 0.0 0.0 
2 8 12.6 14.770 34.000 0.50 -248.7 0.0 0.0 
2 #9 L2.a 14.774 34.000 0.50 -184.7 0.0 0.0 
2 #10 «13.0 14.774 34.000 0.50 -101.3 0.9 0.0 
2m lameel 3.2 14.766 34.000 0.50 .- 4.1 0.0 0.0 
2a 3.14 14.750 34.900 0.50 100.2 0.0 0.0 
2 seek 3/6 14.736 34.000 0.50 100.8 0.0 0.0 
2 #14 «13.7 14.723 34.000 9.50 100.8 0.0 0.0 
2 1S 13.9 14.710 34.000 0.50 100.8 0.0 0.0 
2 16 14.1 14.697 34.000 9.50 100.8 0.0 2.0 
2 WW 46.2 14.684 34.900 0.50 100.8 0.0 0.0 
2 #18 14.64 14.670 34.000 0.50 100.8 0.0 9.0 
2 #19 14.4 14.658 34.000 0.50 100.8 -.00932 -13.497 
2 20) (ee) 14.645 34.000 0.50 100.8 -.01800 -26.005 
Ply = fariaal 14.633 34.000 0.50 100.8 -.02546 -36.398 
2223.8 14.620 34,000 0.50 100.8 -.03118 -43.3855 
2 2238 13.5 14.607 34.000 0.50 100.4 -.03377 -47.921 
2 24 «13.2 14.592 34.000 0.50 100.8 -.03600 -48.452 
. ( 
a] 
“7 TIME (HR) 
© | | | ( ) 
na 
> LATER TEMP. LATER OEPTN 
1%. 3560 5.0473 
I 14.9370 5.3279 
14.9310 27,4124 
S 14.9340 7.8700 
i 14.9420 9.1835 | 
— 14. 9500 4.5379 
4 ‘4.9590 3.8755 
14. 3660 3.1721 ea od 
14.9670 9.4599 Se 
34.3610 3.7736 
a | 14,9490 19.10N2 
oS 14.9270 $0.4379 
: 14.9070 10. 7818 
2 34.8890 11.1138 
ae 1%.8710 11.4199 
14.9540 11.7307 
14.3380 52.022 
14.8213 12.3186 
i=] 14. 3053 12.4846 
=) 14.7300 22.5480 
me 14.7760 12.5274 
ey: +4. 7610 52.4497 
x ta. 74660 12.3462 
= | 4.7269 12.2499 
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Figure 22. Same as figure {l), except 
for Case [Ie with peak at hour 24. 
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OAY UR a r $ TAO QNET woD On 
eoeel 6.1 14.956 34.000 0.50 - 3.5 -.01800 -11.170 
omeee 7.0 14.938 34.000 0.50 -100.8 -.02546 -18.121 
1v3 7.5 14.933 34.000 0.50 -164.3 -.03119 -23.928 = 
14 ee: 14.936 34.000 0.50 +-249.4 -.03477 227.773 & 
1s 3.0 14.944 34.000 0.50 -289.7 - .03600 -29.1249 SS, 
l 6 8.0 14.952 34.000 0.50 -302.5 -.03477 -28.357 —_ = 
1 7 @.1 14.960 34.000 0.50 -289.9 -.03119 -25.656 10 : 
1 3 6.2 14.965 34.000 0.50 -248.7 -.02546 -21.259 
193 a.5 14.964 34.000 0.50 -la4.7 -.01800 -15.397 
1 10 9.8 14.956 34.000 0.50 -101.3 -.00932 - 9.241 
1 iui 9.2 14.940 34.000 0.50 - $.1 0.0 0.0 
1 12 9.5 1$.914 34.000 0.50 100.2 0.0 0.0 
1s 9.9 14.992 34.000 3.50 100.8 0.0 0.0 
1 1é@ 10.2 14.971 34.000 9.50 100.8 0.0 0.0 
1 158 10.5 14.851 34.000 0.50 100.8 0.0 0.0 
1 16 10.8 14.831 34.000 0.50 100.9 0.0 0.0 
eee a Se | 14.9833 34.000 0.50 100.8 0.0 0.0 
119 411.64 14.793 34.0900 0.50 100.8 0.0 0.0 
1419 411.7 14.776 348.000 0.50 100.49 0.0 0.0 
1 20 11.9 14.760 34.000 0.50 100.8 0.0 0.0 a 
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Figure 23. Same as figure (1), except 
for Case IIc with peak at hour S$. 
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